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1  Introduction 

The  purpose  of  the  present  Lecture  Series  is  to  provide  an  up-to-date  review  of  the  experimental  techniques, 
the  theoretical  models,  as  well  as  the  numerical  simulation  strategies  involved  in  the  treatment  of  the  chemical 
character  of  high  temperatures  gases.  To  this  end,  aspects  of  both  materials  and  high  temperatures  uid  sciences 
are  discussed  together,  with  the  aim  of  enhancing  the  improvements  in  the  understanding  of  the  processes  of 
heat  release  on  solid  surfaces  close  to  reactive  gases  taking  place  in  re-usable  hypersonic  ight  vehicles,  bal¬ 
listic  missiles,  and  rockets  exhaust  nozzles.  In  addition,  further  details  are  given  concerning  the  modeling  and 
the  experimental  investigation  of  gas-surface  interactions.  This  latter  aspect  is  strongly  related  to  the  diffusive 
component  of  the  heat  ux  experienced  by  a  vehicle  entering  the  atmosphere.  In  particular,  the  way  in  which 
the  Thermal  Protection  Material  (TPM)  in  uences  the  gas  particle  recombination  on  its  surface  has  a  tremen¬ 
dous  effect  on  the  wall  heat  ux.  As  will  be  shown  later  on,  the  heat  ux  measured  in  the  stagnation  point  of 
a  probe  held  in  a  plasma  wind  tunnel  can  simply  double  when  materials  enhancing  complete  recombination 
are  used  instead  of  inert  material.  Therefore,  the  role  played  by  the  modelling  of  gas/surface  interactions  is  of 
fundamental  importance  in  the  understanding  of  heat-transfer  and  as  a  consequence  in  the  design  of  Thermal 
Protection  Systems  (TPS). 

When  analyzing  the  interaction  between  a  gas  and  a  non  ablating  surface,  at  least  two  approaches  can  be  put  in 
evidence.  A  microscopic  and  a  macroscopic  approach. 

Recent  years  have  seen  a  very  important  progress  in  the  microscopic  theoretical  treatment  of  surfaces  and  pro¬ 
cesses  on  surfaces  described  by  the  so  called  theoretical  surface  science  [31].  The  aim  of  theoretical  surface 
science  is  to  contribute  signi  cantly  to  the  fundamental  understanding  of  the  underlying  principles  governing 
the  geometric  and  electronic  structure  of  surfaces,  together  with  the  processes  taking  place  on  these  surfaces 
such  as  growth  of  surface  layers,  gas-surface  scattering,  friction  or  reaction  at  surfaces  [30].  In  the  framework 
of  atmospheric  (re)entry,  theoretical  surface  science  could  provide  a  very  attractive  tool  for  the  improvement 
of  the  current  techniques  used  to  estimate  the  catalytic  activity  of  thermal  protection  materials.  Some  research 
is  ongoing  in  this  direction  and  published  data  are  already  available  providing  values  of  recombination  prob¬ 
abilities  of  oxygen  on  Silica  based  materials  computed  using  Molecular  Dynamics  techniques  [54].  The  main 
advantage  of  this  microscopic  approach  is  that  a  variety  of  surface  properties  can  be  described  from  rst  prin¬ 
ciples,  i.e.  without  invoking  any  empirical  parameters.  On  the  other  hand,  this  approach,  based  on  expensive 
computations,  has  not  yet  been  used  to  investigate  complicated  o  w  conditions  like  those  typical  of  high  en¬ 
thalpy  facilities.  It  will  certainly  be  worth  to  investigate  this  possibility  in  the  future. 

In  the  eld  of  TPS  design,  a  macroscopic  approach  is  more  widespread,  which  is  closely  coupled  with  the 
use  and  the  simulation  of  high  enthalpy  plasma  o  ws  produced  in  suitable  facilities.  In  this  manuscript,  our 
attention  is  focused  on  this  second  way  of  analyzing  gas-surface  interactions.  In  particular,  the  methodology 
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developed  at  the  Institute  for  Problems  in  Mechanics  (IPM)  has  been  selected  as  a  suitable  instrument  for  the 
determination  of  the  catalytic  activity  of  TPM. 


1.1  IPM  Methodology  to  determine  TPS  materials  catalytic  activity 

Within  the  framework  of  the  IPM  methodology,  the  determination  of  TPS  materials  catalytic  properties  rely 
on  experimental  testing  on  TPS  materials  samples  in  high  enthalpy  facilities  producing  o  w  conditions  close 
to  ight  conditions.  Inductively  Coupled  Plasma  (ICP)  wind  tunnels  are  particularly  suited  for  this  purpose 
because  they  produce  plasma  jets  of  very  high  chemical  purity,  contrary  to  arc  heated  wind  tunnels  in  which 
the  o  w  is  polluted  by  metallic  impurities  originating  from  electrode  erosion,  which  can  contaminate  the  TPS 
material  sample.  This  was  the  motivation  for  the  construction  of  two  ICP  facilities  at  VKI  :  a  large  scale  ICP 
facility  (Plasmatron)  [6]  and  a  small  scale  facility  (Mini torch)  [13]. 

Kolesnikov  et  al.  [39,  34]  proposed  a  methodology  to  estimate  TPS  materials  catalytic  activity.  This  methodol¬ 
ogy  requires  heat  transfer  and  pitot  pressure  measurements  on  one  hand,  and  numerical  simulation  of  the  o  w 
in  the  testing  facility  and  along  the  stagnation  line  in  front  of  the  heat  transfer  probe  on  the  other  hand.  In  the 
following,  we  present  the  main  aspects  of  this  methodology  from  the  analysis  of  which  the  main  motivations  of 
this  contribution  were  inspired. 

The  IPM  methodology  relies  on  a  combination  of  heat  transfer  and  pitot  pressure  measurements  in  a  subsonic 
plasma  jet  produced  in  an  ICP  facility  and  numerical  simulations  of  the  o  w.  In  this  sense  it  represents  a  hybrid 
methodology  which  indeed  builds  on  both  experiments  [15,  27,  40,  12,  2]  and  numerical  simulation  of  high 
enthalpy  o  ws  [75,  43].  Actually,  the  o  w  conditions  in  the  ICP  facility  are  related  to  ight  conditions  using 
the  Local  Heat  Transfer  Simulation  (LHTS)  concept  [35,  3]  so  that  the  measured  heat  ux  es  arc  equivalent  to 
the  actual  stagnation  point  heat  ux  es  of  a  certain  space  vehicle  at  a  point  (altitude/velocity)  of  its  (re)entry 
trajectory.  A  detailed  description  of  the  IPM  methodology  is  presented  in  Ref.  [16]  and  hereafter  we  present 
only  its  main  characteristics. 

This  methodology,  which  is  schematically  represented  in  Fig.  1 ,  is  composed  of  the  following  essential  building 
blocks: 

1 .  An  LTE  with  constant  elemental  fraction  viscous  o  w  simulation  of  the  o  w  in  the  ICP  facility  heating 
chamber  and  in  the  downstream  subsonic  plasma  jet.  The  purpose  of  this  LTE  viscous  o  w  computation 
is  to  compute  the  values  of  the  non-dimensional  parameters  (H,-)  which  characterize  the  o  w  near  the 
stagnation  point  of  the  heat  ux  probe  in  the  plasma  jet. 

2.  Experimental  measurements  of  the  stagnation  point  wall  heat  ux  (qw)  and  of  the  pitot  pressure  (/;)  on  a 
cold  wall  (Tw  ~  300K)  reference  heat  ux  probe  assumed  to  be  fully  catalytic1. 

3.  Reconstruction  of  the  plasma  jet  enthalpy  and  velocity  gradient  at  the  edge  of  the  stagnation  line  using 
the  measured  stagnation  point  heat  ux  and  pitot  pressure. 

4.  Calculation  of  heat  ux  abacus  (optional).  Once  the  stagnation  line  edge  conditions  have  been  deter¬ 
mined,  the  stagnation  line  solver  can  be  run  for  various  wall  catalytic  activity/temperature  combinations 
to  produce  a  heat  ux  abacus  (set  of  qw,  Tw  curves  at  constant  catalytic  activity  yw). 

'For  air  o  ws  the  material  used  is  copper  while  for  COi  mixtures  we  use  Silver. 
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LTE  viscous  flow  simulation 


Experiments 


Set  of  non-dimensional 
— ►  parameters 

n,- 


Pitot  pressure  ( Ap) 
Heat  flux  (qw) 


Reconstruction  of  the  outer  edge  conditions 

Enthalpy  (h§)  Velocity  gradient  (dv/dr|g ) 

|  Non  equilibrium  stagnation  line  computation  ! 
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Outer  edge 
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composition 


—  !  Hydrodynamics  | 
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Computation  of  the  wall  heat  flux  for  various  n  * 
catalytic  efficiency  ( Yw)  /  wall  temperatures 
using  a  non  equilibrium  stagnation  line  solver  0  2 


Y„=* 

-oT~ 

1001 

Yv.-0~~- 

Measurement  on  a 
fully  catalytic 
cold  wall  material 
(Y„=l) 


Test  measurement  on  a  TPM  sample 
""  (Tw,qw) 


Identification  of  Yw 


Figure  1:  Schematic  representation  of  the  IPM  methodology  for  the  determination  of  the  catalytic  activity  of 
TPS  materials 


5.  Experimental  measurement  of  the  stagnation  point  heat  ux  to  a  TPS  material  sample  and  of  its  temper¬ 
ature,  and  determination  of  the  material  catalytic  activity2. 


The  previous  list  shows  clearly  the  highly  hybrid  character  of  the  methodology.  A  particular  attention  should 
be  devoted  to  the  third  step,  usually  referred  as  the  rebuilding  of  the  outer  edge  conditions.  The  logic  on  which 
this  process  is  based  follows  from  the  investigation  of  the  parameters  from  which  the  stagnation  point  heat  ux 
depends.  Such  investigation  reveals  that  these  parameters  are  as  follows  : 


•  the  thermochemical  state  of  the  stagnation  line  outer  edge.  In  the  original  version  of  the  methodology 
the  outer  edge  is  supposed  to  be  under  thermochemical  equilibrium  conditions  de  ned  by  p$,  h§,  and  the 
elemental  fraction  corresponding  to  the  ICP  torch  inlet. 

•  The  radial  velocity  gradient  [(dvldr)s\  at  the  stagnation  line  outer  edge. 

•  The  stagnation  line  nite  thickness  (6)  and  the  product  of  the  normal  (i.e.  axial)  velocity  by  the  normal 
derivative  of  the  radial  velocity  gradient  at  the  stagnation  line  edge  [u$d / dx{dv / dr)§\. 

•  The  wall  temperature  ( Tw )  and  the  catalytic  activity  of  the  sample  surface  with  respect  to  each  of  the 
kth  possible  recombination  processes  (ykw)-  One  of  the  major  assumption  of  the  original  version  of  the 

2This  can  be  done  graphically  using  the  heat  ux  abacus  by  identifying  the  y„,  contour  on  which  the  data  point  (qw,  Tw)  lies,  or 
alternatively  numerically. 
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methodology  consists  in  assuming  all  these  recombination  probabilities  to  be  equal  to  the  same  effective 
probability  yw. 


Among  these  parameters,  only  two  are  unknown  :  the  outer  edge  enthalpy  (I15)  and  the  radial  velocity  gra¬ 
dient  [{dv/dr)s\.  Indeed,  the  stagnation  pressure  (/;„)  and  the  wall  temperature  (Tw)  are  known  experimentally, 
the  recombination  probability  (yjj,  =  yw  ik)  is  equal  to  one  according  to  the  assumption  of  a  fully  catalytic  heat 
ux  probe.  In  addition,  the  stagnation  line  thickness  (d)  and  ugd/dx(dv/dr)s  are  related  to  the  radial  veloc¬ 
ity  gradient  itself  and  to  the  probe  radius  Rlu  through  the  non-dimensional  hydrodynamic  parameters  obtained 
from  the  simulation  of  the  plasma  o  wing  in  the  ICP  facility  performed  assuming  thermochemical  equilibrium 
conditions  with  constant  elemental  fraction  to  be  established  both  in  the  torch  and  in  the  test  chamber.  The  rst 
two  non  dimensional  parameters3  read  as  follows  : 


n 


II2  -  llg 


dx 


(1) 


so  that  the  stagnation  point  heat  ux  can  be  determined  from  the  solution  of  the  stagnation  line  problem  (SL) 
which  implicitly  reads 


n  (dv 
Qwifydi  I 


>  Pdf  III  >  1^2?  Tw). 


(2) 


unknown 


The  additional  equation  needed  to  retrieve  the  two  unknowns  (hg,  (dv/dr)g)  is  given  from  the  knowledge  of  the 
measured  pitot  pressure.  This  measurement  is  related  to  pg  and  ( dv/dr)s  through  a  third  non-dimensional  hy¬ 
drodynamic  parameter  obtained  from  the  LTE  viscous  o  w  simulation  carried  out  assuming  constant  elemental 
fraction,  i.  e. 


n3 


_ Ap _ 

Pd(hs,P6)[Rm(dv/dr)j] 


(3) 


The  two  unknown  o  w  quantities  are  then  determined  by  solving  the  system  (2)-(3)  by  some  suitable  itera¬ 
tive  scheme,  relying  on  the  fact  the  functional  form  (2)  is  provided  implicitly  by  running  a  non-equilibrium 
stagnation  line  solver4. 


In  this  paper  we  wish  to  investigate  and  solve  some  of  the  issues  previously  described  by  merging  and 
exploring  some  of  the  ideas  scattered  in  the  literature,  and  at  the  same  time  present  and  recommend  new 
concepts.  The  outcome  of  this  process  consists  in  the  proposition  of  a  modi  cation  to  the  IPM  methodology  for 
the  determination  of  the  catalytic  properties  of  TPS  materials,  based  on  the  investigation  of  its  main  weak  points. 
The  determination  of  these  weak  points  follows  directly  from  the  description  presented  in  the  previous  section. 
Indeed,  we  notice  that  several  conditions  must  be  satis  ed  in  order  for  the  methodology  to  be  successful : 

3For  the  computation  of  the  non-dimensional  parameters  it  is  necessary  to  de  ne  the  location  of  the  stagnation  line  outer  edge.  This 
position  is  selected  to  be  coincident  with  the  in  ection  point  of  the  axial  velocity  pro  le  along  the  stagnation  line  obtained  from  the 
inductively  coupled  plasma  o  w  simulation.  Further  details  about  the  theoretical  background  on  which  this  choice  is  based  are  given 
in  Ref.  [37],  while  the  low  sensitivity  of  the  II,  to  the  choice  of  8  has  been  veri  ed  in  Ref.  [71] 

4  A  slightly  different  approach,  involving  an  additional  dependent  non-dimensional  parameter,  is  proposed  in  the  original  version  of 
the  IPM  methodology  [16]  and  implemented  in  the  VKI-Rebuilding  code  [28]. 
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•  the  plasma  o  wing  in  the  ICP  facility  should  be  under  thermochemical  equilibrium  and  elemental  diffu¬ 
sion  should  not  affect  the  o  w  sofution. 

•  The  stagnation  line  o  w  should  be  under  chemical  non  equilibrium  to  make  the  wall  heat  ux  sensitive 
to  variations  in  the  surface  catalytic  activity. 

•  The  recombination  probability  of  all  heterogeneous  recombination  processes  should  be  the  same. 


The  investigation  of  the  validity  of  some  of  the  assumptions/conditions  previously  listed  represented  the  main 
motivation  and  objective  of  the  present  study  and  lead  us  to  the  following  main  contributions  of  our  work. 


•  In  the  past,  not  sufficient  credit  has  been  given  to  the  phenomenon  of  elemental  diffusion  within  the 
atmospheric  (re)entry  related  literature.  Indeed,  in  a  wide  variety  of  applications,  this  phenomenon  has 
been  neglected.  In  particular,  this  is  the  case  in  the  1PM  methodology.  The  investigation  of  this  aspect  of 
high  temperature  o  w  modeling  leads  us  to  the  derivation  of  an  explicitly  closed  form  of  the  governing 
equations  of  mixtures  of  reacting  gases  under  LTE  presented  in  Sec.  2.  This  was  possible  reducing 
the  equations  of  chemically  reacting  o  ws  to  an  elegant  system  consisting  of  the  conventional  Navier- 
Stokes  equations  (mass,  momentum,  energy)  complemented  by  an  advection-diffusion  equation  for  the 
mass  fraction  of  each  chemical  element  in  the  mixture.  The  obtained  formalism  is  in  closed  form  in  the 
sense  that  diffusive  ux  es  are  directly  expressed  in  terms  of  gradients  of  the  solution  unknowns,  unlike 
other  formulations  in  which  these  ux  es  are  obtained  in  an  implicit  manner,  by  solving  the  full  system 
of  Stefan-Maxwell  equations.  The  main  contribution  of  this  formalism  consists  in  the  introduction  of 
several  new  LTE  transport  coefficients. 

•  In  Sec.  3  we  present  an  application  of  the  derived  theory  to  the  description  of  CO2/N2  LTE  mixtures 
o  wing  along  a  stagnation  line.  This  leads  to  the  original  computation  of  multicomponent  elemental 
diffusion  coefficients  for  mixtures  containing  more  than  two  eiements. 

•  In  Sec.  4  we  investigate  the  effects  of  efemental  diffusion  on  the  behavior  of  air  and  CO2  plasmas  o  wing 
within  the  VKl-Plasmatron  in  a  certain  range  of  operating  conditions.  In  addition  we  provide  a  charac¬ 
terization  of  hot  air  jet  o  wing  into  the  test  chamber  of  the  VKl-Plasmatron  facility.  This  study  suggest 
a  modi  cation  of  the  original  version  of  the  1PM  methodology  based  on  the  observation  that,  for  suf- 
ciently  high  operating  pressures,  chemical  equilibrium  conditions  with  variable  elemental  fraction  arc 
established  within  the  torch  and  the  test  chamber.  As  a  consequence,  the  stagnation  line  outer  edge  ele¬ 
mental  fractions  should  be  added  to  the  set  of  non-dimensional  parameters  (14,)  needed  in  the  rebuilding 
process. 

•  In  Sec.  5  an  alternative  model  for  the  phenomenological  description  of  gas-surface  interactions  is  pre¬ 
sented  and  tested  in  the  framework  of  TPM  catalycity  determination  for  Mars  entry  applications. 
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2  Theory  of  LTE  viscous  o  ws  with  variable  elemental  fractions 


The  physics  of  collision-dominated  chemically  reacting  o  ws  is  described  by  an  extended  Navier-Stokes  sys¬ 
tem,  consisting  of  the  following  equations  [67]:  global  continuity;  momentum  and  total  energy;  a  separate 
continuity  equation  for  each  species,  including  nite-rate  chemistry;  if  thermal  non  equilibrium  occurs,  an  en¬ 
ergy  equation  for  each  additional  mode  of  freedom  (vibrational,  rotational  and  electronic  energies). 

This  formalism  has  several  drawbacks.  Firstly,  the  equations  are  costly  to  implement  and  solve  numerically. 
Secondly,  many  physical  parameters  essential  for  the  modeling  of  chemistry  and  energy  relaxation  processes 
are  usually  missing.  Thirdly,  even  when  a  numerical  solution  is  obtained,  it  is  by  no  means  obvious  to  interpret 
the  vast  amount  of  information  obtained  (e.g.  concentration  elds  of  a  large  number  of  chemical  species).  For 
these  reasons,  when  chemistry  and  energy  exchanges  are  fast,  it  is  usually  preferable  to  solve  the  more  elegant 
and  less  uncertain  Local  Thermodynamic  Equilibrium  (LTE)  form  of  the  aforementioned  set  of  equations. 

A  major  breakthrough  in  the  eld  of  LTE  o  w  modeling  was  made  by  Butler  and  Brokaw  [9,  7],  who  showed 
that,  assuming  vanishing  diffusive  ux  es  of  chemical  elements,  the  diffusive  transport  of  species  reaction  en¬ 
thalpies  in  the  energy  equation  could  be  incorporated  in  a  straightforward  manner  by  introducing  a  coefficient 
of  ‘thermal  reactive  conductivity’  Ar: 

Nsp 

2  w  A  =  -AVT.  (4) 

S=1 

One  often  makes  use  of  this  result  to  reduce  the  full  set  of  non  equilibrium  equations  to  a  system  formally 
equivalent  to  the  ‘conventional’  Navier-Stokes  equations  (continuity,  momentum  and  energy),  complemented 
by  a  modi  ed  equation  of  state  p{p,  T)  computed  from  statistical  mechanics  assuming  a  x  ed  elemental  com¬ 
position  in  the  o  w.  For  instance,  Vasil’evskii  et  al.  [77]  successfully  used  this  classical  LTE  formalism  to 
simulate  high-pressure  air  inductively  coupled  plasma  o  ws,  imposing  a  21/79  volumetric  ratio  of  oxygen  and 
nitrogen  elements  throughout  the  o  w  eld.  While  appealing  because  of  its  simplicity,  it  is  important  to  under¬ 
stand  that  this  approach  is  approximate  at  best,  since  in  general  the  elemental  composition  varies  signi  cantly 
in  chemically  reacting  o  ws.  Over  the  past  four  decades,  several  LTE  formulations  accounting  for  (de)mixing 
effects  have  been  proposed  in  several  contributions  [48,  23,  24,  25,  26,  55,  70,  32,  22,  72,  73,  78,  68,  41,  45, 
46,  48,  47,  18,  74,  59,  58,  64],  a  description  of  which  is  proposed  in  Ref.  [57].  In  the  remaining  part  of  the 
section  we  will  rstly  present  a  preliminary  analysis  of  the  in  uence  of  elemental  demixing  on  the  behavior 
of  a  mixture  of  reacting  gases  by  computing  the  composition  for  a  x  ed  pressure  and  several  elemental  com¬ 
positions,  in  the  range  of  temperatures  [300  K,  15000  K]  typical  of  applications  to  TPS  testing.  Secondly  we 
will  recall  the  theory  of  LTE  viscous  o  ws  with  variable  elemental  fractions  for  neutral  mixtures  derived  in 
Ref.  [62],  referring  the  extension  to  mixtures  of  ionized  species  to  Ref.  [57].  Although  we  will  focus  on  the 
particular  case  of  air  and  carbon  dioxide  mixtures,  well-suited  for  Earth  and  Mars  entry  applications,  we  wish 
to  point  out  that  the  presented  LTE  formulation  is  applicable  to  any  chemically  reacting  o  w  near  LTE  and  in 
particular  the  hydrocarbon-air  mixtures  used  for  combustion. 


2.1  In  uence  of  elemental  fractions  on  thermochemical  equilibrium  composition 

It  is  well  known  that  the  composition  of  a  mixture  of  reacting  gases  under  thermochemical  equilibrium  can  be 
expressed  as  a  function  of  pressure,  temperature,  and  of  the  elemental  fractions  of  the  elements  shared  among 
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the  mixture  species  [79,  1],  The  contents  of  this  section  deal  mainly  with  the  issue  of  elemental  fraction  vari¬ 
ations  in  reacting  o  ws  under  LTE.  To  attract  the  readers’s  interest,  we  present  a  preliminary  analysis  of  the 
in  uence  of  elemental  fraction  variations  on  the  mixture  composition  for  both  air  and  carbon  dioxide  mixtures. 
As  a  result  we  show  that  a  10%  variation  of  oxygen  elemental  fraction  can  induce  a  50%  difference  in  the 
species  concentration  for  CO2  mixtures. 


T[KJ 


(a)  Mole  fractions  (Z°  =  0.189,  XN  = 
0.811). 


T|K] 


(b)  Mole  fractions  (X°  =  0.21,  XN  = 
0.79). 


T[K] 


(c)  Mole  fractions  (X°  =  0.231,  XN  = 
0.769). 


Figure  2:  In  uence  of  elemental  fractions  on  thermochemical  equilibrium  composition  for  air  and  carbon 
dioxide  mixtures.  (p=l  atm) 

Consider  an  air  mixture  composed  by  the  following  1 1  species  :  N2, 02,  NO,  N,  O,  N} ,  OJ ,  NO+,  N+,  0+,  e". 
Let  us  de  ne  a  reference  mixture  characterized  by  elemental  fractions  X°  =  0.21,  XN  =  0.79,  and  two  addi¬ 
tional  mixtures  obtained  by  perturbing  the  oxygen  fraction  by  an  amount  of  ±10%,  i.e.  X°  =  0.23 1 ,  A'v  =  0.769 
and  X°  =  0.189,  XN  =  0.811.  In  Figs.  2(a),  2(b),  and  2(c)  the  composition  of  these  mixtures  is  shown  as  a 
function  of  temperature  for  a  x  ed  pressure  equal  to  1  atm,  assuming  thermochemical  equilibrium  conditions. 
From  these  gures  appear  that  the  species  mole  fraction  is  little  affected  by  the  elemental  fraction,  whose  in- 
uence  is  basically  limited  to  a  scaling  of  the  various  species  according  to  the  elemental  fraction.  The  reason 
is  because  in  the  air  mixture,  there  are  only  2  species  including  both  N  and  O  atoms,  i.e.  NO  and  NO+,  and  in 
addition  their  concentration  always  remains  small. 


For  carbon  dioxide  mixtures,  the  sensitivity  of  the  chemical  composition  to  elemental  fraction  variations  is 
much  higher  than  for  air.  Indeed,  if  one  considers  an  8  species  mixture  [59]  in  chemical  equilibrium  composed 
by  CO2,  O2,  CO,  C,  O,  C+,  0+,  and  e",  and  computes  the  chemical  composition  as  a  function  of  temperature 
for  a  x  ed  pressure,  very  different  results  are  obtained  depending  on  the  elemental  fraction  used.  Consider  a 
reference  mixture  characterized  by  elemental  fractions  X°  =  |  and  Xc  =  j.  Perturbing  by  ±10%  the  oxygen 
fraction,  two  new  mixtures  are  de  ned:  X°  =  Xc  =  ^  and  X°  =  jjj,  Xc  =  The  mole  fraction  evolution 
as  a  function  of  temperature,  for  a  pressure  of  1  atm,  for  these  three  cases  is  shown  in  Figs.  2(d),  2(e),  and  2(f).  A 
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very  strong  in  uence  of  the  elemental  fraction  on  the  chemical  composition  is  clearly  visible.  Moreover,  from 
the  analysis  of  the  previous  results,  it  follows  that,  in  the  case  of  CO2  mixtures,  the  sensitivity  to  elemental 
fraction  variation  is  much  higher  than  in  the  case  of  air. 

This  will  result  in  an  important  in  uence  of  elemental  demixing  on  heat  ux  in  thermochemical  equilibrium 
CO2  mixtures,  as  will  be  shown  by  the  results  to  be  presented  later  on. 


2.2  Mixtures  of  neutral  species 
2.2.1  Species  ordering  and  nomenclature 


We  represent  mixtures  of  perfect  gases  by  means  of  a  nite  set  of  Nsp  species  S,  amongst  which  we  furthermore 
distinguish  between  Nc  ‘independent  species’  £  consisting  of  pure  elements  (as  such  as  in  their  stable  form) 
and  Nr  ‘combined  species’  R,  for  instance: 


5-species  air  £  =  {O, N},  R  =  { O2 , N2, NO],  and  S  =  R  U  £, 


valid  for  LTE  mixtures  at  pressures  above  0.01  atm  and  for  temperatures  between  300  and  8000  K.  We  accord¬ 
ingly  de  ne  the  three  sets  of  indexes  R  =  { 1, . . . ,  Nr},  £  =  { Nr  +  1, . . . ,  Nsp),  and  S  =  R  U  £.  We  characterize 
the  chemical  composition  of  the  mixture  in  terms  of  mole  fractions  xs  =  njn,  where  ns  and  n  stand  for  the 
molar  densities  of  individual  species  and  of  the  entire  mixture.  Alternatively,  we  can  also  characterize  the  mix¬ 
ture  composition  by  means  of  the  mass  fraction  ys  =  ps/p,  where  ps  and  p  stand  for  the  mass  densities  of  the 
individual  species,  respectively  the  full  mixture. 

We  will  indicate  the  number  of  atoms  of  element  e  contained  in  a  species  s  by  <pes,  for  instance,  for  NO,  0^  =  1 
while  for  O2  (p\  =  2.  This  enables  us  to  de  ne  the  mole  fractions  Xe  and  mass  fractions  Ye  of  elements  in  the 
mixture  as  follows: 


E 

seS 


E  E  <t>esxs  ’ 

eeS  seS 


r 


seS 


Me 

m; 


(5) 


where  Ms  is  the  molar  mass  of  species  5,  related  to  the  mixture  molar  mass  M  =  X  xsMs.  We  introduce  the 
diffusion  velocity  V  v  with  respect  to  the  mass-averaged  velocity  of  the  mixture  u.  The  molar  and  mass  ux  es 
of  species  s  are  then  given  by  respectively  Ws  =  ns\s  and  Jv  =  Msns\s. The  mole  and  mass  ux  es  of  element 
e  are  then  given  by  _ 

Ne  =  £#Wj;  Je  =  Me Ne  (e  e  £).  (6) 

seS 


We  use  bold  type  fonts  to  indicate  vectors  in  the  physical  space.  With  ‘bar  notation’,  we  refer  to  arrays 
containing  species,  elemental  or  reactive  properties,  with  respective  lengths  of  Nsp,  Nc  and  Nr. 
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2.2.2  Computation  of  composition  under  LTE 


Species  continuity  equations 


We  consider  the  commonly  encountered  o  w  regime  in  which  chemical  reactions  are  relatively  rare  w.r.t.  elastic 
collisions,  such  that  they  do  not  play  an  important  role  in  the  thermalization  of  species  in  the  o  w  (unlike  the 
‘kinetic  chemical  equilibrium  regime’  considered  in  the  rst  part  of  Ref.  [18],  for  which  chemical  reactions 
and  elastic  collisions  are  treated  on  the  same  level).  The  concentration  of  each  species  may  then  be  determined 
from  a  respective  species  continuity  equation  [67] : 

dt(pys)  +  V-(py.vu)  +  V-(M/ Ws)  =  u)s  (7) 


where  u  stands  for  the  mass-averaged  velocity  of  the  mixture  and  a>s  is  the  mass  production/destruction 
term  [79]  of  species  s  due  to  chemical  reactions.  The  number  ux  es  of  species  respect  the  mass  conserva¬ 
tion  constraint  _ 

MsWj  =  0  (8) 

seS 

and  obey  the  Stefan-Maxwell  equations  [22,  32] 


M  xsW k  -  xkWs 

~p  h  VMU 


(9) 


for  s  £  S.  To  keep  the  analysis  as  simple  as  possible,  we  have  neglected  effects  of  pressure  and  thermal 
diffusion;  note  however  that  these  could  be  included  without  any  particular  problem  as  shown  in  Ref.  [57]. 


Elemental  continuity  equations  for  neutral  mixtures 


As  pointed  out  by  Murphy,  we  need  to  solve  additional  element  advection-diffusion  equations  to  determine  the 
elemental  composition  of  the  mixture  [46,  58].  The  solution  of  these  equations  will  then  be  used,  together  with 
two  independent  thermodynamics  variables,  temperature  and  pressure  for  example,  to  compute  the  mixture 
composition  as  the  solution  of  a  non-linear  system  [79].  To  obtain  these  equations,  we  multiply  Eqs.  (7)  by 
< pesMe/Ms  and  sum  over  all  species.  Since  no  elements  are  created  in  the  considered  chemical  reactions,  the 
mass  fraction  of  any  element  e  obeys  the  following  equation: 

dt(pYe)  +  V-(puYe)  +  V-Je=0  (eeB).  (10) 

In  Ref.  [62]  we  have  shown  that  under  LTE  conditions,  the  mass  diffusion  ux  of  elements  J e  can  be  expressed 
in  terms  of  gradients  of  elemental  mass  fractions  and  temperature.  For  simplicity,  we  consider  o  ws  at  constant 
pressure,  consistent  with  the  neglect  of  pressure  diffusion  in  Eqs.  (9).  Once  again,  we  remark  that  this  does 
not  imply  any  fundamental  limitation  and  pressure  diffusion  could  be  included  easily  if  needed  (see  Ref.  [57]). 
Following  the  approach  proposed  in  Ref.  [62],  the  elemental  diffusive  ux  es  are  expressed  introducing  the 
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elemental  multicomponent  diffusion  coefficients  Def  and  thermal  demixing  coefficients  DTe.  This  lead  to  the 
following  expression  for  the  elemental  diffusive  ux  es  : 

Je  =  ~Yj  PDeFYf  -  pD^VT.  (11) 

/e£ 

The  second  term  in  the  right  hand  side  of  the  above  result  will  in  general  generate  nonzero  elemental  diffusion 
ux  es  even  when  the  initial  elemental  composition  is  uniform. 


2.2.3  Diffusive  transport  of  enthalpy 


In  reacting  o  ws,  the  diffusion  of  species  affects  the  mixture  energy  balance  through  the  heat  ux  term: 

qrf  =  ^WA-  (12) 

seS 

In  a  general  non  equilibrium  case,  to  compute  qj,  one  should  determine  all  the  Wv  as  a  solution  of  Eq.  (9)  and 
then  compute  the  above  linear  combination,  as  done  in  the  methodology  of  Refs.  [70], [18]  and  [58].  On  the 
other  hand,  under  thermochemical  equilibrium,  Eq.  (12)  can  be  cast  under  a  particular  form  which  avoids  the 
computation  of  the  Ws.  A  rst  step  in  this  direction  was  made  by  Butler  and  Brokaw  [9],  who  showed  that, 
under  the  assumption  of  vanishing  elemental  ux  es,  the  diffusive  heat  ux  takes  the  form  of  Eq.  (4).  In  this 
section,  we  extend  the  work  of  Butler  and  Brokaw  to  the  more  general  case  when  elemental  ux  es  are  nonzero. 
As  shown  in  Ref.  [62],  for  a  mixture  of  reacting  gases  under  thermochemical  equilibrium,  at  constant  pressure, 
the  diffusive  heat  ux  is  proportional  to  both  temperature  and  elemental  concentration  gradients,  and  can  be 
expressed  as: 

2  WA  =  -Ur  +  A)VT  -  ^  4  Lvr.  (13) 

sE»S  eE& 

From  Eq.  (13),  we  see  that  in  an  LTE  o  w,  the  diffusive  heat  ux  vector  consists  of  three  different  parts: 


(1)  the  ‘thermal  reactive  conductivity’  coefficient  Ar  is  identical  to  the  well-known  results  by  Butler  and 
Brokaw  [9,  7]  and  takes  into  account  diffusive  transfer  of  species  enthalpies  in  the  absence  of  elemental 
diffusion. 

(2)  The  ‘thermal  demixing  conductivity’  coefficient  Ad  corrects  for  the  additional  diffusive  heat  transfer  that 
occurs  due  to  nonzero  elemental  diffusive  ux  es  when  elemental  mass  fraction  gradients  are  zero. 

(3)  Finally,  the  ‘elemental  heat  transfer  coefficients’  AeEL,  take  into  account  heat  transfer  due  to  elemental 
demixing  driven  by  gradients  in  elemental  composition. 


In  the  next  section  we  will  present  an  application  of  this  formulation  to  investigate  diffusion  phenomena  along 
the  stagnation  line  present  in  front  of  a  probe  placed  in  the  test  chamber  of  an  inductively  coupled  plasma 
facility. 
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3  Diffusion  phenomena  in  reacting  mixtures  :  applications  to  stagnation  line 
o  ws 


This  section  is  inspired  by  the  research  presented  in  Ref.  [63],  where  the  authors  present  a  detailed  analysis  of 
diffusion  phenomena  in  o  ws  of  reacting  mixtures  along  a  stagnation  line  present  in  front  of  a  probe  placed 
in  the  test  chamber  of  an  inductively  coupled  plasma  facility.  In  the  following  we  limit  our  description  to  the 
presentation  of  the  results  of  the  solution  of  the  stagnation  line  differential  problem  and  we  refer  the  interested 
reader  to  our  previous  publications  [59,  58,  57]  for  a  detailed  presentation  of  the  stagnation  line  equations.  In 
these  references  we  also  discuss  the  models  used  for  the  computation  of  the  transport  and  thermodynamics 
properties  of  reacting  mixtures  such  as  those  used  to  obtain  the  following  results. 

The  following  mixture  has  been  selected  as  representative  for  Mars  entry  : 

Mars-8  species  :  £  =  [C,  N,  O),  %  =  {C02,  CO,  02,N2,  NO). 


The  computations  carried  out  using  the  Mars-8  mixture  represent  a  direct  application  of  the  theory  presented  in 
Sec.  2.  This  represents  the  rst  attempt  to  simulate  mixtures  of  hetero-nuclear  molecules  with  more  than  two 
elements  using  the  closed  form  of  the  equations  derived  in  Ref.  [62].  The  use  of  the  theory  presented  in  Sec.  2, 
will  moreover  allow  to  show  the  improvements  brought  by  the  theory  itself  in  terms  of  physical  understanding 
of  the  o  w  behavior.  Although  stagnation  line  o  ws  are  often  under  chemical  non  equilibrium  conditions,  our 
previous  analysis  [58]  showed  that  for  a  CO2  mixture,  assuming  a  fully  catalytic  wall  provides  an  estimation  of 
wall  heat  ux  close  to  the  one  obtained  by  assuming  LTE  conditions  to  be  established  along  the  stagnation  line, 
provided  that  elemental  demixing  is  taken  into  account.  In  addition,  the  elemental  fraction  pro  les  obtained 
under  chemical  non  equilibrium  follow  closely  those  under  chemical  equilibrium.  Thus  we  will  present  the 
solution  of  the  LTE  stagnation  line  problem  formulated  following  the  equations  presented  in  Ref.  [57]. 

To  present  our  analysis  of  the  stagnation  line  o  w  we  will  start  de  ning  the  test  case  to  be  investigated  and  then 
the  solution  of  the  stagnation  line  o  w  will  be  presented  with  strong  emphasis  on  the  description  of  diffusion 
phenomena. 


3.1  Test  case  de  nition 

The  o  w  conditions  considered,  characteristic  of  Martian  entry  [66],  are  presented  in  Table  1. 

The  presented  outer  edge  conditions  and  geometrical  parameters,  have  been  determined  from  a  detailed  numer¬ 
ical  [33]  study  of  the  o  w  eld  inside  the  inductively  coupled  plasma  generator  and  test  chamber  of  the  VKI’s 
plasmatron  wind  tunnel  [5].  This  allows  for  the  computation  of  some  non-dimensional  parameters  related  to 
the  stagnation  line  edge  in  terms  of  the  boundary  layer  thickness,  the  axial  velocity,  and  the  velocity  gradient 
and  further  details  about  their  de  nition  are  available  in  previous  publications  [16]. 

In  the  following  a  single  chemical  regime  is  considered,  corresponding  to  LTE  with  variable  elemental  fraction 
(LTE-VEF),  where  the  o  w  is  in  thermochemical  equilibrium  and  the  composition  is  computed  as  a  function 
of  pressure,  temperature,  and  local  elemental  fractions,  obtained  from  the  solution  of  the  set  of  elemental  con¬ 
tinuity  equations.  Two  formulations  have  been  used. 
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Table  1 :  Operating  Conditions. 


Tw  [K] 

300 

TS\K\ 

5827 

Ps  [Pa] 

7000 

R,n  [mm] 

25 

5  [mm] 

9 

vsd/dy(dug/dx)/(dug/dx)2  [-] 

1.77 

dug/dx  [s-1] 

913 

LTE-VEF  (1)  we  obtain  the  elemental  concentration  by  solving  Eqs.  (10),  where  the  element  diffusive  ux  es 
are  computed  as  a  linear  combination  of  the  species  ones  obtained  as  the  solution  of  the  Stefan-Maxwell 
equations  [Eq.  (9)].  The  computed  species  diffusive  ux  es  are  then  used  in  the  energy  equation  to  evaluate 
the  diffusive  transport  of  enthalpy  [58]  (X  Jshs)- 

LTE-VEF  (2)  we  obtain  the  elemental  concentration  by  solving  Eqs.  (10),  where  the  thermal  demixing  and 
multicomponent  diffusion  coefficients  are  used  to  compute  the  element  diffusive  ux  es  [Eq.  (1 1)].  More¬ 
over,  an  alternative  form  of  the  energy  balance,  based  on  the  use  of  Eq.  (13)  is  used,  introducing  the 
correction  to  Ar  as  well  as  the  elemental  heat  transfer  coefficients. 


3.2  Stagnation  line  solution 

In  this  section  we  present  the  solution  of  the  stagnation  line  equations  for  the  conditions  speci  ed  in  Tab.  1 .  We 
start  our  analysis  by  discussing  the  temperature  and  enthalpy  pro  les  shown  in  Fig.  3(a). 

There  we  present  the  results  obtained  for  the  LTE-VEF  regime  using  the  two  formulations  LTE-VEF  (1)  and 
LTE-VEF  (2)  previously  de  ned.  As  observed,  the  results  obtained  with  the  two  formulations  are  identical. 
The  same  match  is  observed  for  all  o  w  variables  along  the  stagnation  line,  supporting  the  correctness  of  the 
alternative  formulation  recalled  in  Sec.  2. 

In  Fig.  3(b)  we  present  the  species  concentration  pro  les.  Starting  from  CO,  O,  C,  and  N  at  the  outer  edge,  we 
see  how  their  concentration  decreases  providing  a  mixture  of  CO2,  O2,  N2  at  the  wall  and  a  small  amount  of 
NO.  From  the  analysis  of  this  picture,  it  appears  that  elemental  demixing  tends  to  reduce  the  amount  of  carbon 
and  nitrogen  at  the  wall  while  it  enhances  the  oxygen  concentration.  This  is  re  ected  by  the  presence  of  an  ex¬ 
cess  of  O2  in  the  low  temperature  region,  which  would  have  been  zero  if  the  wall  elemental  fractions  had  been 
the  same  as  at  the  outer  edge.  This  behavior  is  con  rmed  by  the  pro  les  presented  in  Fig.  3(c),  where  we  plot 
the  difference  between  the  local  elemental  mass  fractions  and  their  outer  edge  values  {Ye(y)-Ye\g\.  In  addition, 
by  looking  to  the  element  concentration  pro  les  [Fig.  3(c)]  we  notice  that  the  lack  of  carbon  and  nitrogen  at  the 
wall  is  compensated  by  a  higher  concentration  around  y/6  =  0.4,  while  for  oxygen  we  observe  a  minimum  in 
this  position  and  a  subsequent  concentration  rise  as  the  wall  is  approached. 

Going  further  with  the  analysis  of  the  pro  les  presented  in  Fig.  3(c),  we  notice  that  they  present  a  non-monotone 
behavior  and  hereunder  we  give  an  alternative  explanation  for  this  observation.  To  discuss  the  elemental  mass 
fraction  behavior  we  focus  on  the  elemental  continuity  equations  [Eq.  (10)].  It  is  clear  how  the  knowledge  of 
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(c)  Difference  between  the  elemental  mass  fractions  and 
their  outer  edge  values  (Ye(y)-Ye\6). 


Figure  3:  Comparison  of  the  computed  results  along  the  stagnation  line  obtained  using  the  LTE-VEF  (1)  and 
the  LTE-VEF  (2)  formulations. 


Je  will  help  understanding  the  evolution  of  Ye  knowing  that  the  velocity  is  aligned  with  the  stagnation  line  and 
directed  towards  the  wall. 

For  this  purpose  we  present  the  pro  les  of  J~e  in  Fig.  4  for  the  three  elements  contained  in  the  mixture.  As 
expected  [Eq.  (10)]  we  observe  a  correspondence  between  the  local  extrema  of  J~e  and  those  of  Ye.  In  addition, 
we  notice  that  the  diffusive  ux  es  of  elements  are  zero  all  over  the  rst  10%  of  the  stagnation  line  and  they 
start  to  grow  around  y/6  ~  0.1.  This  shows  that  over  the  rst  10%  of  the  stagnation  line,  Eq.  (10)  simpli  es  to 
Je  =  0,  meaning  that  convection  is  by  far  negligible.  Moving  further  towards  the  outer  edge  a  local  extremum 
is  encountered,  followed  by  an  increase  for  C  and  N,  while  for  oxygen  a  decrease  of  JO  is  observed.  From  this 
we  notice  that  in  the  region  de  ned  by  y/6  >  0.1,  the  convective  term  will  differ  from  zero  and  will  counteract 
the  derivative  of  the  pro  les  presented  in  Fig.  4.  At  this  point,  it  is  important  to  observe  that,  although  the 
pro  les  of  Fig.  4  allow  the  interpretation  of  the  elemental  fraction  pro  les  presented  in  Fig.  3(c),  it  is  by  no 
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Figure  4:  Elemental  diffusive  ux  es  (-^Te). 


means  evident  to  give  an  explanation  for  the  behavior  of  the  'Je  if  they  are  computed  as  a  linear  combination  of 
the  solutions  of  Eqs.  (9)  as  done  in  the  formulation  LTE-VEF  (1)  [58,  64].  Now,  to  improve  our  understanding, 
we  use  the  theory  introduced  in  Sec.  2  and  applied  in  the  LTE-VEF  (2)  formulation  of  the  equations. 

Indeed,  to  understand  the  behavior  of  the  element  diffusive  ux  es  we  discuss  now  the  two  components  of  Je. 
The  rst  proportional  to  the  temperature  gradient  (pDjVT)  and  the  second  obtained  as  a  linear  combination 
of  the  element  mass  fractions  gradients  (p£  DeqVYq).  In  Fig.  5(a),  we  present  the  rst  part  of  the  element 
diffusive  ux  es.  From  Fig.  3(a)  we  notice  that  dT  j dy  >  0  V y  and  therefore  the  sign  of  pD^VT  will  depend  only 
on  Dj,  p  being  obviously  a  positive  quantity. 

As  shown  in  Ref.  [62,  63]  and  con  rmed  in  Fig.  5(b),  the  DTe  change  their  sign  as  temperature  rises,  for  x  ed 
pressure  and  elemental  fractions.  In  the  present  case,  the  pressure  is  constant  along  the  stagnation  line  but  the 
elemental  fraction  varies.  Therefore,  both  temperature  and  elemental  composition  will  in  uence  the  transport 
properties  and  therefore  concur  to  the  establishment  of  the  pro  les  presented  in  Fig.  5(a). 

Close  to  the  wall,  the  thermal  demixing  coefficients  are  zero  and  do  not  contribute  to  the  element  diffusive 
ux  es.  As  temperature  rises,  we  observe  an  increase  in  module  of  all  contributions  until  a  maximum  is  reached. 
Then  we  notice  a  decrease  and  a  subsequent  change  in  sign  typical  of  the  thermal  demixing  coefficients  behav¬ 
ior. 

Indeed,  from  Fig.  5(b),  we  see  how  the  behavior  of  the  thermal  demixing  coefficients  is  similar  to  the  one  of 
the  rst  contribution  to  the  element  diffusive  ux  es. 

In  Fig.  5(c),  we  plot  the  second  contribution  to  the  elemental  diffusive  ux  es.  Also  this  contribution  starts 
from  zero  at  the  wall  and  follows  a  non  monotone  behavior  passing  through  a  change  in  sign  for  O  and  C.  The 
change  in  sign  observed  for  the  second  contribution  to  the  elemental  diffusive  ux  es  is  not  due  to  the  elemental 
multicomponent  diffusion  coefficients  but  to  the  elemental  mass  fractions.  Indeed,  as  shown  in  Ref.  [63]  and 
observed  in  Fig.  5(d),  the  multicomponent  diffusion  coefficients  have  the  same  sign  in  the  temperature  range  of 
interest.  This  allows  for  the  following  interpretation  of  diffusion  phenomena  along  the  stagnation  line.  Since 
temperature  decreases  from  the  outer  edge  towards  the  wall,  the  temperature  gradient  is  always  positive  along 
the  stagnation  line.  In  sufficiently  low  temperature  regions,  where  chemical  reactions  do  not  occur,  the  thermal 
demixing  coefficients  are  zero  and  do  not  induce  demixing.  As  a  consequence,  since  elemental  multicomponent 
diffusion  coefficients  are  not  zero  at  low  temperatures,  the  elemental  composition  stays  constant  ensuring  mass 
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(c)  Second  contribution  to  the  element  diffusive  uxes  (d)  Multicomponent  elemental  diffusion  coefficients. 

(pE^vi”). 


Figure  5:  Elemental  diffusive  ux  components,  elemental  multicomponent  and  thermal  diffusion  coefficients. 

conservation  (^Te|w  =  0).  As  we  move  away  from  the  wall  (y  <  0.05),  the  thermal  demixing  coefficients  causes 
a  non  zero  contribution  to  the  elemental  diffusion  ux  which  tends  to  decrease  oxygen  concentration  ( D TQ  >0), 
increase  nitrogen  concentration  (DTN  <  0),  and  slightly  decrease  but  then  increase  again  carbon  concentration 
(DTC  >  0  and  DJC  <  0).  As  soon  as  variations  in  elemental  composition  are  induced  (VFe  +  0),  a  counteracting 
contribution  to  the  element  diffusive  ux  es  appear  which  tends  to  smooth  out  elemental  concentration  pro  les 
until  the  outer  edge  is  reached.  In  the  central  part  of  the  stagnation  line,  the  two  contributions  to  element 
ux  es  interact  until  the  point  in  which  the  thermal  demixing  coefficients  change  sign.  There,  to  counteract  this 
demixing  term,  the  sign  of  element  concentration  gradients  needs  to  change  since  the  multicomponent  diffusion 
coefficients  have  a  constant  sign  [Fig.  5(d)].  The  counteracting  character  of  the  term  £  pDeqWYq  with  respect  to 
pD^VT  justify  its  label  of  ‘mixing  term’  and  the  comparison  of  Figs.  5(a)  and  5(c)  clearly  shows  this  behavior. 
The  mathematical  origin  of  this  diffusive  character  of  the  mixing  term  has  to  be  searched  in  the  non  negativity 
of  the  matrix  containing  the  elemental  multicomponent  diffusion  coefficients.  To  prove  the  latter  property  of  the 
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matrix  Deq  we  computed  its  eigenvalues  along  the  stagnation  line.  Since  this  matrix  is  singular  [ Deq  =  0], 
one  eigenvalue  is  k\  =0  and  the  remaining  two  (k^,kj)  have  been  found  to  be  positive  and  monotonically 
increasing  as  a  function  of  y  [63].  This  result  shows  the  non  negative  character  of  the  matrix  Deq. 

As  a  summary,  we  wish  to  point  out  that,  from  the  knowledge  of  the  set  of  transport  properties  including  both 
elemental  thermal  demixing  and  multicomponent  diffusion  coefficients,  we  can  justify  the  evolution  of  the 
diffusive  ux  es  and  therefore  understand  the  shape  of  the  elemental  fractions  pro  les  one  can  observe  in  the 
solution  of  the  stagnation  line  problem. 


Figure  6:  Fleat  ux  components. 

We  now  move  to  the  analysis  of  the  heat  and  mass  transfer  along  the  stagnation  line  by  means  of  conduction 
and  diffusion.  As  shown  in  Sec.  2,  under  the  assumption  of  local  thermodynamic  equilibrium,  the  diffusive 
transport  of  enthalpy  can  be  split  into  two  parts  related  to  VT  and  VYq  respectively.  As  a  consequence,  as 
shown  in  Eq.  (2),  the  computation  of  the  heat  ux  involves  the  determination  of  two  terms  one  proportional  to 
the  temperature  gradient  [-(d  +  Ar  +  Ap)VT  ]  and  the  other  proportional  to  the  linear  combination  of  element 
concentration  gradients  [-  £  AqE[yYq).  Moreover,  the  heat  ux  in  a  mixture  of  reacting  gases  is  composed  of  a 
conductive  part  (qc  =  -AVT),  to  which  the  diffusive  transport  of  enthalpy  needs  to  be  added  (qj  =  £  hsJs).  To 
investigate  the  heat  ux  along  the  stagnation  line  and  highlight  the  relative  importance  of  the  various  contribu¬ 
tions  we  start  presenting  the  evolution  of  q,  qc,  and  q,/  along  the  stagnation  line  in  Fig.  6.  There  we  observe  that 
the  total  heat  ux  presents  an  increasing  monotone  behavior  starting  from  the  outer  edge  until  the  wall,  acting 
to  balance  the  convection  of  enthalpy.  At  the  same  time  we  notice  that  this  smooth  evolution  appears  to  be  the 
result  of  a  quite  complex  shape  of  the  two  contributions  qc  and  q d  presented  in  the  same  gure.  This  highly 
non-linear  behavior  of  the  two  heat  ux  components  is  the  result  of  the  evolution  of  six  transport  coefficients 
(/l,  Ar,  Ad,  AqEI)  and  four  gradients  (V7\  VYq).  To  provide  an  exhaustive  description  of  the  heat  ux  along  the 
stagnation  line,  we  start  investigating  the  behavior  of  the  three  thermal  (d),  thermal  reactive  (Ar)  and  demixing 
(Ad)  conductivities  shown  in  Fig.  7(a). 

As  expected,  we  notice  Ar  to  be  the  major  contribution  to  the  rst  part  of  the  heat  ux  followed  by  A.  We  also 
notice  that  Ad  has  a  non-monotone  behavior  including  changes  in  sign,  highlighting  the  fact  that  demixing  tends 
to  increase  or  decrease  heat  ux  depending  on  the  local  temperature  and  elemental  fractions  as  already  observed 
in  Ref.  [64].  Moreover,  in  terms  of  heat  ux  intensity,  we  notice  that  around  y/6  =  0.3,  Ad  ~  A  showing  how, 
neglecting  Ad  with  respect  to  A  could  lead  to  important  errors.  The  knowledge  of  the  evolution  of  the  transport 
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(a)  Thermal,  reactive,  and  demixing  conductivities  along 
the  stagnation  line. 


(b)  Heat  ux  components  (AKVT,  Ar,VT)  along  the  stag¬ 
nation  line. 


10s 


Figure  7:  Heat  ux  components,  thermal  conductivities,  and  elemental  heat  transfer  coefficients. 


coefficients  presented  in  Fig.  7(a)  helps  us  understanding  the  behavior  of  the  various  heat  ux  contributions. 
Indeed,  the  non-linear  behavior  of  qc  can  be  easily  explained.  A  is  an  increasing  monotone  function  from  the 
wall  towards  the  outer  edge,  while  dT / dy  starting  from  the  outer  edge  increases,  decreases  and  then  increases 
again  as  shown  in  Fig.  3(a).  Therefore  qc  follows  mainly  the  evolution  of  dT /dy.  For  what  concerns  q d,  we 
need  to  consider  separately  the  ve  contributions  -ArVT,  —ApVT,  and  -AqF[VYq. 

In  Fig.  7(b),  we  present  the  evolution  of  the  rst  two  terms  and  their  sum.  The  evolution  of  these  two  contri¬ 
butions  is  more  complex  than  the  previous  one  since  the  two  associated  transport  coefficients  present  a  highly 
non-linear  behavior  with  extrema  in  different  positions  than  the  temperature  gradient.  As  far  as  the  elemental 
heat  transfer  coefficients  are  concerned,  we  present  their  evolution  in  Fig.  7(c).  There  we  notice  that  they  all 
have  a  non-monotone  behavior  re  ecting  the  highly  reacting  character  of  the  o  w.  These  transport  coefficients 
lead  to  three  contributions  to  the  heat  ux  presented  in  Fig.  7(d),  which  are  of  the  same  order  of  magnitude  for 
the  three  elements.  The  sum  of  these  contributions  is  also  plotted  in  the  same  gure  to  show  the  importance  of 
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this  term  with  respect  to  the  total  heat  ux. 

From  the  analysis  of  the  pro  les  presented  in  Figs.  7(b)-7(d),  we  can  justify  the  evolution  of  concurring  to 
the  determination  of  the  nal  shape  of  the  heat  ux.  As  a  nal  step  we  analyze  in  detail  the  relative  importance 
of  the  various  terms  contributing  to  the  heat  ux.  In  Fig.  8,  we  present  the  results  of  this  analysis  showing  the 
evolution  of  the  four  following  ratios  along  the  stagnation  line  :  -( AdT /dy)/q,  -(ArcIT /dy)/q,  -{Aj^dT / dy)/q , 
and  -(2  AqELdYq / dy) / q. 


Figure  8:  Normalized  heat  ux  components. 

The  analysis  of  these  pro  les  help  us  to  quantify  the  importance  of  each  term.  Indeed,  this  shows  that  the 
major  contribution  to  the  heat  ux  comes  from  the  thermal  reactive  conductivity  (4^),  the  second  comes  from 
conduction  (A)  and  the  two  remaining  contributions  are  of  the  same  order  of  magnitude.  The  analysis  of  the 
curve  -{ApdT / dy) / q  reveals  that  the  contribution  related  to  elemental  demixing  is  higher  than  the  one  due  to 
thermal  conductivity  for  0.1  <  y/6  <  0.3.  The  two  contributions  -A[>dT / dy  and  -  £  AqELdYq /dy  present  almost 
an  opposite  behavior  for  y/6  <  0.4.  On  the  other  hand,  for  y/6  >  0.4  we  notice  that  -  £  AqELdYq/dy  becomes 
higher  than  -A^dT /dy  approaching  the  10%  of  the  total  heat  ux  around  y/6  ~  0.6  and  decreasing  to  5%  at  the 
outer  edge.  The  almost  opposite  behavior  of  the  two  latter  contributions  clearly  shows  that  elemental  diffusive 
ux  es  have  not  a  very  important  in  uence  on  the  total  heat  ux.  Indeed  if  AJe  =  0  Ve  e  £,  the  Butler  and 
Brokaw  thermal  reactive  conductivity  is  enough  to  compute  qf/.  For  the  conditions  previously  analyzed,  the 
presence  of  non  vanishing  elemental  ux  es  acts  on  the  mixture  energetic  behavior  with  a  heat  ux  contribution 
which  is  at  most  of  the  order  of  5%  of  the  total  heat  ux. 


4  Applications  to  ICP  o  w  computation 


In  this  section  we  present  a  detailed  analysis  of  chemical  equilibrium  and  non  equilibrium  plasma  o  ws  both  in 
the  torch  and  in  the  test  chamber.  In  the  remaining  part  of  the  section  we  limit  our  discussion  to  the  presentation 
of  the  results  of  several  plasma  o  w  simulations  and  we  refer  the  interested  reader  to  previous  publications  [64, 
75,  17]  for  the  description  of  the  ICP  model  used  to  obtain  these  results.  In  the  same  references  the  reader 
will  nd  an  extensive  description  of  the  models  used  to  compute  the  thermodynamic  and  transport  properties 
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of  reacting  mixtures  such  as  those  used  to  obtain  the  results  to  be  presented  shortly. 

In  Sec.  4.1  we  will  present  results  concerning  only  the  torch  domain  for  carbon  dioxide  plasmas,  while  in 
Sec.  4.2  we  extend  our  investigation  to  the  hot  air  jet  o  wing  into  the  chamber. 

The  two  following  mixtures  will  be  considered  in  this  section  to  simulate  air  and  carbon  dioxide  plasmas  : 

Air-11  species  :  &  =  {0,N,e+},K  =  {02,  N2,  NO,  NO+,  0+,  N+,  0\,  N+},  and  S  =  k  U  £, 

C02-8  species  :  6  =  {0,C,e+},R  =  {C02,  02,  CO,  0+,  C+),  and  S  =  K  U  £. 

The  analysis  presented  in  the  rst  part  of  this  section  is  based  on  the  study  presented  by  Rini  et  al.  in  Ref.  [65] 
where  a  code-to-code  validation  is  presented  for  carbon  dioxide  plasmas  in  the  VKI  plasmatron  torch.  To  this 
end  the  authors  present  the  results  of  the  simulation  of  CO2  plasma  o  ws  within  the  VKI  plasmatron  torch 
carried  out  obtained  independently  from  the  two  groups  with  the  VKI-ICP  and  the  IPM-cr  codes.  The  results  of 
this  comparison  showed  how  the  1  -D  approach  of  the  a  code  fails  in  the  description  of  the  Electric  eld  within 
the  torch.  In  addition,  the  importance  of  elemental  diffusion  on  the  behavior  of  CO2  plasma  o  ws  is  assessed 
and  only  this  part  will  be  presented  in  the  following  whose  motivation  is  recalled  hereafter: 


(1)  accurate  carbon  dioxide  plasma  o  w  simulations  are  rare  in  the  literature  related  to  aerospace  sciences. 
Moreover,  up  to  the  best  author  knowledge,  no  simulation  is  available  where  elemental  demixing  is 
considered  under  thermochemical  equilibrium  conditions.  For  this  reason  we  investigate  the  phenomena 
of  elemental  diffusion  in  LTE  carbon  dioxide  plasmas  to  evaluate  its  effects  on  the  o  w  behavior. 


The  second  section  is  devoted  to  the  physico-chemical  characterization  of  the  jet  o  wing  in  the  test  chamber 
for  air  plasmas,  whose  motivation  is  twofold  : 


(1)  the  methodology  proposed  by  IPM  [35]  to  estimate  the  catalytic  properties  of  thermal  protection  system 
materials  supposes  thermochemical  equilibrium  conditions  to  be  established  in  the  torch  and  in  the  test 
chamber.  However,  there  is  no  proof  that  equilibrium  conditions  exist  and  a  detailed  analysis  assessing  the 
quality  of  this  assumption  is  not  available.  To  supply  this  need  we  present  both  chemical  non  equilibrium 
and  chemical  equilibrium  computations  to  provide  a  range  of  pressure  in  which  equilibrium  computations 
are  as  accurate  as  those  obtained  under  chemical  non  equilibrium,  using  two  nite  rate  chemistry  models. 

(2)  As  a  result,  we  will  be  able  to  assess  the  in  uence  of  the  nite  rate  chemistry  model  and  operating 
pressure  on  the  o  w  behavior  both  in  the  torch  and  especially  in  the  test  chamber. 

4.1  Analysis  of  carbon  dioxide  plasma  o  ws  within  the  torch 

We  now  move  to  the  analysis  of  the  numerical  results  obtained  for  an  operating  condition  of  the  VKI  plasmatron 
representative  of  Mars  entry  [33].  The  analysis  presented  in  this  section  is  inspired  on  the  study  conducted  by 
Rini  et  al.  in  Ref  [65] ,  omitting  the  part  concerning  the  code-to-code  comparison.  Therefore,  in  the  following, 
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we  only  present  the  investigation  of  elemental  diffusion  to  assess  the  importance  of  the  variations  in  carbon  and 
oxygen  elemental  concentrations  on  the  o  w  behavior.  In  Tab.  2,  we  present  the  operating  condition  selected 
as  representative  of  Mars  entry  [33]  used  to  produce  the  results  to  be  presented  shortly. 


Table  2:  Operating  conditions  for  the  VKI-Plasmatron  operating  with  carbon  dioxide. 


Ambient  and  wall  temperature  [K]  300 

Power  injected  into  the  plasma  [kW]  90 

Frequency  [MHz]  0.37 

Inlet  Swirl  [deg]  0 

Mass  o  w  [g/s]  8 

Operating  pressure  [Pa]  7000 


The  o  w  has  been  simulated  under  LTE  with  both  constant  and  variable  elemental  fractions  and  both  com¬ 
putations  have  been  earned  out  on  the  same  mesh  of  102  by  92  cells  for  which  grid  convergence  is  assured 
(Ref.  [64]).  All  results  presented  in  the  following  are  converged  at  least  10  order  of  magnitude  (based  on  the 
drop  of  the  L?  residual  norm). 

To  investigate  the  in  uence  of  the  elemental  fraction  variation  on  the  o  w  behavior,  we  used  the  upgraded  ver¬ 
sion  of  the  VKI-ICP  code  [57]  to  simulate  both  the  LTE-VEF  and  the  LTE  with  Constant  Elemental  Fraction 
(LTE-CEF)  formalisms.  For  that  purpose,  the  plasma  in  the  torch  has  been  considered  both  under  LTE-VEF 
and  LTE-CEF  conditions.  For  the  operating  conditions  presented  in  Tab.  2  we  obtain  the  elemental  distribution 
presented  in  Fig.  9(a)  in  terms  of  carbon  elemental  mole  fraction  (Xc). 

In  Fig.  9(b),  we  present  the  temperature  contours  and  in  Fig.  9(c)  the  o  w  pattern,  both  within  the  torch,  ob¬ 
tained  under  LTE-VEF  conditions.  In  these  gures,  we  observe  the  presence  of  a  quite  wide  high-temperature 
region  where  the  maximum  temperature  is  around  ~  10  000  K,  as  expected  for  the  sped  ed  power  injected  in 
the  o  w.  The  mixture  enters  the  torch  with  a  reference  elemental  composition  of  2/3- 1/3  for  oxygen  and  carbon 
respectively.  Entering  the  torch,  the  mixture  is  subject  to  very  strong  temperature  gradients  which  produce  a 
non  zero  demixing  contribution  (-pDjVT)  to  the  elemental  diffusion  ux  [62],  inducing  elemental  composi¬ 
tion  variations.  Once  a  non-uniform  elemental  composition  eld  is  established,  a  counteracting  mixing  term 
(YjpDckVYk),  proportional  to  the  elemental  composition  gradients,  tends  to  re-establish  a  uniform  elemental 
eld.  These  two  mechanisms  act  together  with  convection  to  provide  the  nal  steady  elemental  composition 
pattern.  From  the  analysis  of  Fig.  9(a),  at  least  three  zones  can  be  observed.  The  rst  one  represents  an  almost 
cylindrical  zone  delimited  by  7.5  <  r  <  8.0  cm  close  to  the  inlet  and  6.0  <  r  <  8.0  cm  at  the  outlet.  There, 
we  do  not  observe  substantial  variations  in  the  elemental  fraction  and  this  is  mainly  due  to  the  fact  that,  in  that 
region,  the  o  w  is  continuously  lied  by  a  stream  coming  from  the  inlet  with  the  reference  elemental  composi¬ 
tion.  In  addition,  the  wall  temperature  being  x  ed,  the  axial  temperature  gradient  is  negligible  close  to  the  wall 
leading  to  the  absence  of  diffusion  in  the  same  direction. 

A  second  zone  is  visible  around  the  torch  axis  both  inside  and  after  the  recirculation  bubble.  There,  we  ob¬ 
serve  that  carbon  has  the  tendency  to  diffuse  out  from  hot  regions  towards  lower  temperatures,  leading  to  an 
axial  zone  poorer  in  carbon  than  in  the  case  where  a  constant  elemental  composition  is  arbitrarily  imposed. 
A  third  region  is  visible  at  the  intermediate  radii  and  a  bit  downstream  of  the  inlet.  There,  the  carbon  mole 
fraction  reaches  the  maximum  value  of  0.38,  i.e.  ~  13%  more  than  the  inlet  value.  This  zone  of  high  carbon 
concentration  is  therefore  the  result  of  the  diffusion  of  carbon  from  the  high  temperature  region  around  the  axis 
accompanied  by  the  difficulty  for  carbon  to  diffuse  towards  the  quartz  tube  where  the  concentration  is  mainly 
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(b)  Temperature  (K). 
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(c)  Stream  function. 

Figure  9:  Carbon  dioxide  plasma  o  ws  within  the  torch  (LTE-VEF  regime). 


x  ed  by  the  cold  mixture  convected  from  the  inlet  to  the  torch  exit. 

The  previous  description  of  the  carbon  elemental  eld  is  re  ected  by  the  pro  le  depicted  in  Fig.  10(a)  which 
corresponds  to  the  torch  exit.  There  we  clearly  see  that  the  carbon  pro  le  has  a  non-monotone  behavior.  As 
we  approach  the  exit,  the  o  w  eld  became  almost  axial  and  radial  convection  can  be  neglected  as  shown  in 
Fig.  9(c).  This  behavior  of  the  outlet  carbon  elemental  fraction,  observed  for  nitrogen  in  air  mixtures  [64],  could 
be  predicted  by  observing  the  evolution  of  the  elemental  diffusion  coefficients  for  the  carbon  element  [59]  as  a 
function  of  temperature.  In  Figs.  10(b)  and  10(c),  we  present  the  in  uence  of  elemental  fraction  variations  on 
the  mixture  composition  and  on  the  mixture  temperature.  From  the  analysis  of  Fig.  10(b),  a  big  in  uence  of 
elemental  fraction  variations  on  the  species  concentration  is  observed.  Starting  from  the  wall  we  observe  that, 
since  the  carbon  molar  fraction  is  slightly  smaller  than  1/3,  a  small  amount  of  O2  is  present  close  to  the  wall 
as  a  result  of  a  LTE-VEF  computation.  Moving  further  towards  the  axis  the  differences  between  the  CEF  and 
the  VEF  computations  become  more  evident.  Indeed,  the  decrease  of  Xc  up  to  a  local  minimum  around  6  cm 
causes  a  slower  CO2  dissociation  accompanied  by  the  displacement  to  lower  radii  of  both  the  O2  and  the  CO 
peaks.  Then,  because  of  the  successive  CO  dissociation  we  observe  the  formation  of  both  atomic  carbon  and 
oxygen,  where  the  second  one  is  more  affected  by  the  elemental  fraction  variation  because  of  the  retardation  of 
both  O2  and  CO  peaks.  Finally,  atomic  species  decrease  their  concentration  close  to  the  axis  because  of  ioniza¬ 
tion  due  to  the  high  temperature  around  10000  K.  On  the  axis  we  therefore  nd  a  carbon  fraction  of  ~  0.275, 
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(a)  Carbon  elemental  mole  fraction. 


Figure  10:  Radial  pro  les  at  the  outlet  of  the  torch. 


i.e.  ~  16%  less  than  the  reference  value. 

The  nal  in  uence  of  elemental  demixing  is  observed  in  Fig.  10(c),  where  the  outlet  temperature  pro  les  for 
both  CEF  and  VEF  are  compared.  We  notice  the  temperature  to  be  overestimated  in  the  case  of  constant  ele¬ 
mental  composition  starting  from  r  <  7  cm.  By  comparing  Fig.  10(c)  with  the  results  presented  in  Ref.  [64],  we 
easily  realize  that  the  in  uence  of  elemental  diffusion  on  the  temperature  pro  le  is  much  higher  for  CO2  than 
for  air  mixtures. 


4.2  Characterization  of  air  plasma  o  w  in  the  test  chamber 

We  now  move  to  the  analysis  of  air  plasma  o  ws.  The  results  to  be  presented  shortly  concern  both  the  in¬ 
vestigation  of  the  torch  and  the  hot  exhaust  jet  in  the  chamber.  As  discussed  at  the  beginning  of  this  section, 
the  forthcoming  results  will  allow  to  verify  for  which  operating  conditions  of  the  VKI  plasmatron,  the  nu¬ 
merical  simulation  of  LTE  plasma  o  ws  is  as  accurate  as  the  chemical  non  equilibrium  one,  both  in  the  torch 
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and  especially  in  the  test  chamber.  The  consequences  of  this  nding  will  have  a  direct  impact  on  the  range  of 
applicability  of  the  methodology  used  for  the  estimation  of  catalytic  properties  of  thermal  protection  system 
materials  [35], 

The  mesh  used  for  the  solution  of  the  o  w  eld  equations  is  divided  into  two  blocks:  the  rst  one  (60  by  44 
cells)  is  used  to  discretize  the  domain  within  the  torch,  while  the  second  block  (60  by  67  cells)  reproduces  a 
part  of  the  test  chamber.  All  results  presented  in  the  following  are  converged  at  least  10  order  of  magnitude 
(based  on  the  drop  in  the  Li  residual  norm). 

The  operating  conditions,  selected  as  representative  of  Earth  (re)entry  [51],  are  recalled  in  Tab.  3. 


Table  3:  Plasmatron  operating  conditions. 


Ambient  and  wall  temperature  [K] 

350 

Power  injected  into  the  plasma  [kW] 

80 

Frequency  [MHz] 

0.37 

Inlet  Swirl  [deg] 

0 

Mass  o  w  [g/s] 

8 

Operating  pressures  [kPa] 

5,  10,  20,  30 

For  each  simulation,  the  three  chemical  regimes  chemical  non  equilibrium  (CNEQ),  LTE-CEF,  and  LTE-VEF 
have  been  considered. 

In  the  following  we  present  a  qualitative  analysis  based  on  the  investigation  of  the  contour  plots  of  the  main 
o  w  properties  and  we  refer  the  interested  reader  to  Ref.  [57]  for  a  more  detailed  description  of  the  o  w  behav¬ 
ior  based  on  the  analysis  of  temperature  and  elemental  molar  fraction  pro  les  extracted  in  several  crucial  parts 
of  the  computational  domain. 

We  rstly  address  the  issue  of  the  in  uence  of  nite  rate  chemistry  on  the  results  of  the  simulation.  To  this 
end,  two  models  have  been  selected  as  representative  of  Earth's  atmosphere  chemistry  (Park  [53]  and  Dunn- 
Kang  [29]).  Several  numerical  experiments  have  been  conducted  to  assess  how  the  operating  pressure  drives 
the  o  w  behavior  when  one  of  the  two  models  is  used.  As  a  result,  we  notice  that  the  in  uence  of  the  nite  rate 
chemistry  is  only  limited  to  low  pressures  (around  5000  Pa).  Indeed,  as  the  operating  pressure  is  set  to  higher 
values,  the  results  obtained  under  CNEQ  with  both  models  are  in  very  good  agreement. 

We  start  presenting  temperature  and  elements  concentration  both  in  the  torch  and  in  the  chamber  by  means  of 
contour  plots.  Two  of  the  four  operating  pressures  speci  ed  in  Tab.  3  have  been  selected  to  perform  this  study, 
the  lowest  (5000  Pa)  and  the  highest  (30000  Pa). 

Focusing  on  these  two  pressures,  we  will  be  able  to  easily  show  how  the  differences  observed  in  the  o  w  solu¬ 
tion,  due  to  differences  in  the  reaction  rates,  are  visible  only  in  the  low  pressure  case  and  disappear  when  the 
operating  pressure  is  increased  up  to  30000  Pa.  The  results  to  be  presented  shortly  have  been  obtained  under 
LTE-VEF  conditions  and  under  CNEQ  using  two  nite  rate  chemistry  models:  the  Park  model  [53],  indicated 
as  CNEQ-P,  and  the  Dunn- Kang  model  [29]  named  as  CNEQ-D&K. 

Each  non  equilibrium  result  is  compared  with  the  respective  one  obtained  with  the  other  nite  rate  model  but 
also  with  the  result  of  the  simulation  carried  out  under  thermochemical  equilibrium  with  variable  elemental 
fraction.  This,  at  the  same  time,  allows  for  the  assessment  of  the  in  uence  of  the  nite  rate  chemistry  model 
and  to  answer  the  question  of  whether  or  not  equilibrium  computations  give  the  same  results  as  non  equilibrium 
simulations  for  high  pressures  in  the  chamber. 

In  Figs.  11(a)  and  11(b)  we  plot  the  temperature,  respectively  the  oxygen  elemental  mass  fraction  both  in  the 
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(a)  Temperature  (K)  contours  obtained  with  two  nit  e 
rate  models. 


(b)  Oxygen  elemental  molar  fraction  obtained  with  two 
nite  rate  models. 


CNEQ-P 


(c)  Comparison  of  the  temperature  (K)  contours  ob¬ 
tained  with  the  Dunn-Kang  model  (upper)  and  under 
LTE-VEF  (lower). 


(d)  Comparison  of  the  temperature  (K)  contours  ob¬ 
tained  with  the  Park  model  (upper)  and  under  LTE-VEF 
(lower). 


(e)  Comparison  of  the  oxygen  elemental  molar  fraction 
contours  obtained  with  the  Dunn-Kang  model  (upper) 
and  under  LTE-VEF  (lower). 


(f)  Comparison  of  the  oxygen  elemental  molar  fraction 
contours  obtained  with  the  Park  model  (upper)  and  un¬ 
der  LTE-VEF  (lower). 


Figure  1 1 :  Air  plasma  o  w  within  the  VKI  plasmatron  facility  under  CNEQ,  LTE-CEF,  and  LTE-VEF  condi¬ 
tions  (5000  Pa). 


torch  and  in  the  chamber.  The  analysis  of  these  gures  reveals  the  presence  of  two  zones  in  which  the  compar¬ 
ison  of  the  results  bring  us  to  different  conclusions.  Indeed,  within  the  torch,  we  notice  that  the  CNEQ-P  and 
CNEQ-D&P  predictions  differ  considerably  leading  to  different  temperatures  and,  to  a  lesser  extent,  elemental 
fractions.  From  Fig.  11(a)  we  observe  that  the  Dunn-Kang  model  predicts  a  slightly  wider  hot  (T>  6000  K) 
region  in  the  torch  acting  on  the  structure  of  the  plasma  ball.  As  soon  as  we  move  towards  the  torch  exit,  we 
start  noticing  a  rapid  dwindling  of  the  disparity  between  the  two  predictions  which,  at  least  from  a  qualitative 
point  of  view,  show  the  same  temperature  eld  in  the  jet,  corresponding  to  the  most  interesting  part  of  the 
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(a)  Temperature  (K)  contours  obtained  with  two  nit  e 
rate  models. 


(b)  Oxygen  elemental  molar  fraction  obtained  with  two 
nite  rate  models. 
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(c)  Comparison  of  the  temperature  (K)  contours  ob¬ 
tained  with  the  Dunn-Kang  model  (upper)  and  under 
LTE-VEF  (lower). 


(d)  Comparison  of  the  temperature  (K)  contours  ob¬ 
tained  with  the  Park  model  (upper)  and  under  LTE-VEF 
(lower). 


CNEQ-P 


LTE-VEF 


(e)  Comparison  of  the  oxygen  elemental  molar  fraction 
obtained  with  the  Dunn-Kang  model  (upper)  and  under 
LTE-VEF  (lower). 


(f)  Comparison  of  the  oxygen  elemental  molar  fraction 
obtained  with  the  Park  model  (upper)  and  under  LTE- 
VEF  (lower). 


Figure  12:  Air  plasma  o  w  within  the  VKI  plasmatron  facility  under  CNEQ,  LTE-CEF,  and  LTE-VEF  condi¬ 
tions  (30000  Pa). 


o  w  eld  from  the  point  of  view  of  TPM  testing. 

For  what  concerns  the  oxygen  elemental  fraction,  we  notice  that  the  nite  rate  chemistry  model  has  little  in- 
uence  on  the  elements  concentration  eld,  showing  once  again  that  demixing  is  not  a  non  equilibrium  effect 
[Fig.  1 1(b)].  In  Figs.  1 1(c)- 1 1(d)  and  1  l(e)-l  1(f),  we  present  a  comparison  between  the  results  obtained  under 
chemical  non  equilibrium  with  the  Park  and  Dunn-Kang  models  with  those  obtained  under  chemical  equilib¬ 
rium  with  variable  elemental  fractions.  In  Figs.  1 1  (c)- 1 1  (d)  we  focus  on  the  temperature  contours  while  in 
Figs.  1  l(e)-l  1(f),  we  investigate  the  oxygen  elemental  fraction.  As  expected  [64],  the  temperature  contours  ob- 
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tained  under  chemical  non  equilibrium  differ  considerably  from  those  obtained  with  the  equilibrium  formalism 
and  this  for  both  nite  rate  chemistry  models.  From  the  oxygen  elemental  fraction  presented  in  Figs.  11(e)- 
1 1(f)  we  notice  a  certain  difference  in  the  prediction  of  the  elemental  concentration  eld  between  equilibrium 
and  non  equilibrium  even  though  the  order  of  magnitude  of  the  oxygen  variations  is  the  same. 

We  now  move  to  the  analysis  of  the  higher  pressure  (30000  Pa)  and  we  start  our  analysis  presenting,  as  done 
for  the  lower  pressure  case,  the  comparison  between  non  equilibrium  results.  In  Figs.  12(a)-12(b)  we  present 
the  temperature,  respectively  the  elemental  fraction  contours  both  in  the  torch  and  in  the  chamber.  From  the 
analysis  of  the  temperature  contours,  we  notice  that  at  an  operating  pressure  of  30000  Pa,  the  results  obtained 
with  the  two  models  arc,  at  least  from  a  qualitative  point  of  view,  equivalent.  The  same  conclusion  is  drawn 
after  inspection  of  Fig.  12(b)  where  attention  is  focused  on  the  oxygen  elemental  fraction.  Differences  between 
the  prediction  with  the  Park  and  the  Dunn-Kang  models  are  hardly  visible.  In  the  following,  we  present  the 
comparison  between  the  previous  results  with  the  respective  equilibrium  computations.  The  critical  analysis 
of  these  comparison  will  enhance  the  importance  of  the  previous  conclusions.  In  Figs.  12(c)-12(d)  and  12(e)- 
12(f)  we  present  the  direct  comparison  between  the  temperature,  respectively  the  elemental  mole  fraction  con¬ 
tours  obtained  both  under  chemical  non  equilibrium  and  equilibrium  with  variable  elemental  fraction.  From 
Figs.  12(c)-12(d)  we  realize  that  the  temperature  contours  predicted  by  the  simulation  of  equilibrium  condi¬ 
tions  coincide  with  those  obtained  under  chemical  non  equilibrium  for  both  nite  rate  chemistry  models  after 
the  last  coil.  Indeed,  some  small  differences  are  visible  within  the  torch  until  the  last  coil  is  reached.  To  further 
decrease  these  differences  in  the  torch  we  just  need  to  increase  the  operating  pressure.  On  the  other  hand,  at 
the  torch  exit,  but  most  importantly  in  the  test  chamber,  a  perfect  match  is  observed.  Similar  conclusions  can 
be  drawn  for  the  elemental  composition  after  inspection  of  Figs.  12(e)- 12(f). 

The  results  just  described  have  a  direct  impact  on  the  methodology  used  for  the  estimation  of  the  catalytic  prop¬ 
erties  of  thermal  protection  materials.  Indeed,  thanks  to  the  previous  analysis  we  can  argue  that  for  sufficiently 
high  pressures  the  equilibrium  formalism  allows  for  the  determination  of  results,  both  in  the  torch  and  in  the 
chamber,  characterized  by  an  accuracy  very  close  to  the  one  we  can  achieve  with  non  equilibrium  simulations. 
This  conclusion  will  allow  researchers  performing  inductively  coupled  plasma  o  w  simulations  for  the  charac¬ 
terization  of  TPM,  to  obtain  accurate  results  using  an  equilibrium  formulation  at  a  fraction  of  the  cost  required 
to  simulate  the  same  conditions  under  non  equilibrium. 


5  Practical  application  to  TPM  testing 


The  von  Karman  Institute  has  been  involved  since  1997  in  the  determination  of  catalytic  properties  of  ther¬ 
mal  protection-system  materials.  Recently,  with  the  initiation  of  several  Mars  exploration  projects,  interest  has 
been  attracted  to  Mars  entry  problems  [19,  20,  11].  The  two  inductively  coupled  plasma  (ICP)  facilities  of 
the  Institute  have  been  equipped  to  work  with  CO2  mixtures  [59].  The  VKI  high-enthalpy  o  w  solvers  arc 
interfaced  with  libraries  to  compute  thermodynamic,  transport,  and  non  equilibrium  chemistry  properties.  The 
thermodynamic  and  transport  libraries  have  been  recently  updated  to  incorporate  CO2  mixtures.  In  this  section 
the  physico-chemical  description  of  a  CO2  o  w  is  completed,  allowing  for  simulations  in  chemical  non  equilib¬ 
rium,  characterized  by  both  bulk  and  surface  reactions.  A  critical  output  of  the  LHTS  methodology  is  the  heat 
load  on  the  stagnation  point  of  a  ying  body,  where  large  gradients  of  temperature  and  mass  concentrations  are 
present.  These  features  make  this  point  a  suitable  benchmark  for  testing  thermochemical  models.  The  o  w  eld 
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in  the  vicinity  of  the  stagnation  point  has  been  analyzed  by  means  of  stagnation-line  o  w  calculations.  These 
computations  are  applied  for  the  determination  of  surface  catalytic  efficiency  in  CO2  mixtures,  leading  to  the 
determination  of  so-called  heat-  ux  maps.  The  analysis  that  ensues  is  based  on  the  study  presented  by  Rini  et 
al.  in  Ref.  [59]  whose  motivation  is  recalled  here  after  : 


(1)  the  simulation  of  gas/surface  interactions  in  CO2  o  ws  requires  a  special  care  because  of  the  coupling 
between  wall  reactions  sharing  oxygen  atoms.  This  is  correctly  described  by  the  model  presented  by  Rini 
et  al.  in  Ref.  [59],  which  will  be  successfully  applied  to  describe  heterogeneous  recombination  in  carbon 
dioxide  mixtures. 

(2)  The  model  proposed  in  Ref.  [59]  for  the  description  of  gas/surface  interactions  is  compared  against  the 
one  presented  by  Scott  [69]  showing  that  the  latter  one  is  not  consistent  with  element  mass  conservation. 

(3)  In  the  aerospace  related  literature,  the  usual  de  nition  of  a  ‘catalytic  surface’  is  often  given  as  the  prop¬ 
erty  of  a  certain  material  to  promote  recombination  of  atoms  or  molecules  on  the  surface  itself.  As  a 
consequence,  a  fully  catalytic  surface  is  said  to  be  the  one  that  forces  a  complete  recombination  induc¬ 
ing  local  conditions  which  tend  to  the  equilibrium  limit.  Following  this  observation,  several  researchers 
working  in  the  eld  of  TPS  design  treat  the  gas/surface  interactions  for  fully  catalytic  materials  by  com¬ 
puting  the  local  mixture  composition  as  the  one  corresponding  to  equilibrium  conditions  de  ned  by  the 
local  pressure  and  temperature  but  using  a  reference  elemental  fraction,  usually  corresponding  to  the  free 
stream.  Thanks  to  the  analysis  of  the  results  to  be  presented  shortly,  we  will  show  that  to  de  ne  the 
equilibrium  conditions  representing  the  limit  toward  which  a  fully  catalytic  condition  should  tend,  the 
local  elemental  fraction  should  be  used  instead  of  the  free  stream  one. 

(4)  A  crucial  aspect  of  the  methodology  for  the  estimation  of  the  catalytic  properties  of  TPM  is  the  heat  ux 
map  expressing  the  relation  between  wall  heat  ux,  wall  temperature  and  the  effective  recombination 
probability.  As  a  nal  summary  of  the  previous  investigation,  we  will  show  how  the  outer  edge  elemental 
fraction  affects  the  heat  ux  map.  This  will  allow  for  the  estimation  of  the  in  uence  that  elemental 
demixing  in  the  jet  could  have  on  the  prediction  of  the  effective  recombination  probability  of  thermal 
protection  materials. 


5.1  Stagnation  line  o  w  computations 

In  this  section,  we  present  the  solution  of  the  stagnation  line  differential  problem  for  a  carbon  dioxide  mixture 
under  chemical  non  equilibrium  conditions.  As  in  Sec.  3,  our  attention  will  be  focused  mainly  on  the  presen¬ 
tation  of  the  results  while  we  refer  the  interested  reader  to  our  previous  publications  [59,  4]  for  a  complete 
description  of  the  stagnation  line  equations.  The  test  case  de  nition  is  based  on  several  simulations  and  experi¬ 
ments  presented  by  Kolesnikov  in  Refs.  [39]  and  [38]  referring  to  the  100-kW  IPG-4  plasmatron  (IPM).  Many 
calculations  have  been  carried  out  at  VKI  [56]  covering  a  wide  range  of  enthalpy  conditions.  To  investigate 
the  importance  of  bulk  and  surface  chemistry,  we  focus  our  attention  on  two  operating  conditions  characterized 
by  different  outer  edge  enthalpies  presented  in  Tab.  4,  where  li$  is  the  mixture  enthalpy  at  the  outer  edge,  Vs 
the  incoming  o  w  velocity  and  p  the  operating  pressure.  This  will  lead  to  different  chemical  compositions 
and  temperatures  at  the  outer  edge.  As  a  consequence,  different  chemical  processes  will  be  enhanced  by  the 
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bulk  recombination  and  the  presence  of  atomic  carbon  at  the  outer  edge  in  the  second  case,  will  in  uence  the 
de  nition  of  the  set  of  wall  reactions. 


Table  4:  Flow  conditions. _ 

Low  enthalpy  case  High  enthalpy  case 


hs  [MJ/kg] 

15.3 

39 

Ps  [atm] 

5.8-KT2 

0.1 

Rm  [mm] 

25 

25 

8  [mm] 

10 

10 

v6d/dy(du6/dx)/(du6/dx)2  [-] 

2.09 

2.09 

dug/dx  [s 1] 

1794 

2558 

The  o  w  is  considered  to  be  under  equilibrium  conditions  at  the  boundary  layer  edge.  Therefore,  the  outer  edge 
temperature  is  the  result  of  an  iterative  calculation  in  which  the  thermodynamic  modeling  of  the  o  w  plays  an 
important  role.  The  corresponding  values  of  temperature  are  5922  K  and  8127  K  for  the  low  and  high  enthalpy 
cases  respectively. 

Table  5:  Chemical  reactions  used  in  the  Stagnation  Line  computations.  (M=C02,  CO,  O2,  O,  C) 

Dissociation  reactions 
C02  +  M  CO  +  O  +  M 
CO  +  M^C  +  O  +  M 
02  +  M^±0  +  0  +  M 
Neutral  exchange  reactions 

co+o^c+o2 

co2+o^co+o2 


For  the  given  o  w  conditions  at  the  boundary  layer  edge,  the  heat  ux  at  the  stagnation  point  of  the  model  has 
been  computed  as  a  function  of  the  surface  temperature  Tw  and  the  surface  effective  catalycity  yw,  from  the 
solution  of  the  stagnation  line  equations.  The  analysis  of  the  non  equilibrium  dissociated  subsonic  o  w  around 
the  model  is  based  on  the  assumptions  used  to  derive  the  stagnation  line  equations  [4,  59]  and  the  following 
ones  : 

•  the  gas  is  represented  by  a  5  species  (C02,C0,02,0,C)  viscous  and  heat  conductive  mixture. 

•  The  reaction  set  describing  the  nite  rate  bulk  chemistry  of  the  o  w  is  the  one  described  in  Tab.  5. 

•  The  surface  catalytic  reactions  proceed  according  to  the  Eley-Rideal  mechanism  [21]  and  the  following 
models  are  considered  for  the  various  cases  analyzed. 

2- surface  reactions  model  : 

(1)  O  +  (O  -  S)  ->  02  +  5  (2)  CO  +  (O  -  5)  ->  C02  +  S 

3- surface  reactions  model  : 

(1)  O  +  (O  -  S)  -»  02  +  5  (2)  CO  +  (O  -  S)  -»  C02  +5  (3)  C  +  2(0  -  5)  — »  C02  +  2S 
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Another  important  assumption  in  the  present  model  is  that  all  the  recombining  species  have  the  same  effective 
recombination  probability,  i.e.  for  the  mentioned  wall  chemistry  model  we  have 

rcow  =  yow  =  yew  =  yw-  (14) 

As  the  result  of  the  computations  performed  along  the  stagnation  line,  for  a  50  mm  diameter  probe,  several 
charts  have  been  determined  for  the  stagnation  point  heat  ux  as  a  function  of  Tw  and  yw  e  [0,1],  and  will  be 
presented  later  in  this  section. 


5.1.1  Gas  surface  interactions 


The  description  of  non  equilibrium  o  ws  requires  the  modeling  of  recombination  phenomena  happening  at 
the  solid  surface.  To  this  end  we  have  speci  ed  two  sets  of  heterogeneous  reactions.  The  next  step  consists 
in  modeling  the  way  in  which  the  species  involved  in  these  reactions  exchange  atoms  in  order  to  express  the 
mass  production  term  oj^cat  as  a  function  of  known  quantities.  Scott  [69]  proposes  a  model  of  the  boundary 
conditions  for  the  species  equations  [Eq.  (7)]  based  on  the  following  expression  for  the  net  mass  ux  of  species 
i  at  the  surface 

Nr  Nr  Nsp 

OJi,cat  =  MiM\  ^  vriyr  -  p, j,YM jMlj ,  (15) 

r=l  r= 1  j= 1 


where  yr  is  the  reaction  probability  of  the  r,h  reaction,  M,  the  molar  mass  of  the  ith  species,  and  At.  is  the 
thermal  agitation  ux  of  /-particles  towards  the  surface.  The  matrix  vri  indicates  which  ith  reactants  collide 
upon  the  surface  in  the  rth  reaction,  while  the  matrices  p,/(  de  ne  the  jth  incident  reactants  producing  the  ith 
products  for  the  rth  reaction.  Within  the  methodology  followed  for  the  determination  of  catalytic  properties 
of  TPS  materials,  the  reaction  probability  yr  is  supposed  to  be  independent  from  the  reaction  and  equal  to  an 
effective  recombination  probability  yw. 

Therefore,  the  expression  for  the  impinging  ux  M\,  if  the  Chapman-Enskog  perturbation  term  is  considered, 
reads 


M\  = 


■  7w 


kTw 
2  nM: 


This  formulation  ensures  a  zero  net  total  mass  ux  at  the  wall  but  does  not  satisfy  in  general  the  element 
conservation,  indeed,  considering  a  catalytic  surface  model  constituted  by  the  two  reactions  O  +  O  — >  O2  and 
CO  +  O  — »  CO2,  one  notices  a  coupling  between  the  two  chemical  processes.  This  in  some  way  establishes  a 
link  between  the  reaction  probabilities  of  the  two  reactions  which  is  not  speci  ed  in  the  previous  formulation. 
Elaborating  further  on  this  consideration,  an  alternative  approach  is  proposed.  Let  us  de  ne  a  rate  of  reaction 
per  unit  surface  and  unit  time  for  the  wall  reaction^,  which  represents  the  frequency  at  which  reactions  take 
place  on  the  surface. 

With  this  in  mind,  a  formalism  similar  to  the  one  used  for  bulk  reactions  is  followed.  The  number  of  moles 
produced  per  unit  surface  and  unit  time  are  given  by 


a>j 


Mi 


(16) 
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where,  for  v.  +  0 


M; 


Xr  =  7; 


v. 


(17) 


The  surface  reaction  rate  Xr  is  therefore  related  to  the  recombination  probability  of  species  i  involved  in  the  rth 
reaction  y'  by  the  relation  (17).  As  the  reaction  rate  Xr  has  to  be  independent  of  the  species  i,  it  results  that 
the  species  reaction  probabilities  y[  are  linked  by  the  constraints  y-Aif  / v'jr  —  constant  =  Xr,  for  all  species  such 
that  v'ir  +  0. 

In  this  case,  the  ux  of  particle  Afj"  is  given  by  Aij  =  2/(2 -y/)n,-  V 'kTw/(2nMj ),  while  y,-  =  £r=i,ivr  yj  represents 
the  recombination  probability  of  species  i.  The  computation  of  the  boundary  condition  for  the  species  continuity 
equations  will  therefore  consists  in  rstly  calculate  Xr  and  then  evaluating  the  production  term  by  means  of 
Eq.  (16).  The  puipose  of  this  alternative  formulation  is  not  to  better  describe  the  physical  phenomenon  of 
heterogeneous  catalycity,  which  is  hidden  in  the  de  nition  of  the  different  recombination  probabilities.  But  at 
least,  it  provides  a  model  which  is  consistent  with  both  the  global  mass  and  element  conservation  at  the  wall, 
as  will  be  demonstrated  by  numerical  experiments.  Moreover,  it  should  be  noticed  that  the  quantity  Xr  depends 
only  on  the  reaction  considered  and  not  on  the  species  involved  in  it,  which  represents  another  similitude  with 
the  bulk  reaction  mechanism. 


5.1.2  Low  enthalpy  case 


We  start  our  analysis  considering  the  low  enthalpy  case  for  which  only  the  2-surface  reactions  model  has  been 
applied,  using  both  Scott’s  model  [69]  and  the  formulation  presented  in  the  previous  section.  This  reaction 
scheme  is  based  on  the  recombination  of  atomic  oxygen  and  carbon  oxide  to  produce  molecular  oxygen,  re¬ 
spectively  carbon  dioxide.  The  success  of  these  two  processes  is  based  on  the  simultaneous  presence  of  O  and 
CO  close  to  the  surface.  The  atomic  oxygen  will  then  concur  to  both  recombination  processes  while  CO  only 
to  one  of  them.  The  description  of  the  kinetics  of  these  two  reactions  requires  the  knowledge  of  the  rates  of 
each  chemical  process.  These  frequencies  will  be  computed  on  the  basis  of  the  effective  probability  that  an 
atom  of  oxygen  has  to  be  involved  in  the  ith  process  (y'0)  and  on  the  probability  that  a  molecule  of  CO  has  to 
be  involved  in  the  j,h  reaction  (yJco)-  Obviously  these  effective  recombination  probability  are  related  by  mass 
conservation  and  constrained  by  Eq.  (14).  Nevertheless,  a  correct  treatment  of  these  issues  is  missed  in  several 
gas/surface  interaction  models  [69,  35],  whereas  they  are  correctly  described  by  the  formalism  presented  in 
Sec.  5.1.1  which  in  the  present  case  reads  : 

To  +  To  =  7o,  7co  =  7co,  720Ml0  =  y2coAilco,  (18) 

where,  thanks  to  Eq.  (14),  yo  =  yco  =  TV  Eqs.  (18)  constitute  a  3  X  3  system  in  the  unknowns  ylQ,  y2Q,  y2CQ 
which  must  also  satisfy  the  inequalities  o  <y\  <  yw.  The  solution  is  easily  obtained 


7Co  =  7  iv, 


2  2 
7o  =  7co 


^ CO 

Mi’ 


7o  =  7h>  -  To- 


(19) 


This  solution  is  valid  as  long  as  y2Q  <  yw.  If  it  is  not,  then  it  means  that  one  of  the  assumptions  made  does 
not  hold  (e.g.  equal  recombination  probability  for  all  species).  In  practice,  y2Q  was  always  below  yw  for  our 
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computations.  Then,  the  reaction  frequencies  [Eq.  (17)]  are  obtained  as 


Xi  = 


y'oK 


Xi 


-  TcoM 


l 


2  ’  r coy,'co > 

leading  to  the  following  expressions  for  the  species  production  terms 


co2 


=  y2coH 


CO-  Wo  =  ~y2ooMco 


o)o2 


yloK 


’  M, 


02 


toe 

Wc 


=  o, 


wo 

M0 


(20) 


=  “(To  +  ya>K-  (21) 


(a)  McKenzie  bulk  chemistry  model.  (b)  Park  bulk  chemistry  model. 

Figure  13:  Species  mass  fraction  and  non  dimensional  temperature  pro  les  [Surface  reaction  model  presented 
in  Eq.  (21),  /t,5=15.3MJ/kg,  r„,=300K],  [0:  wall;  1:  5], 

We  now  move  to  the  presentation  of  the  solution  of  the  stagnation  line  problem  [59].  All  the  computations 
showed  in  the  following  are  converged  for  a  100  points  discretization  of  the  stagnation  line.  A  grid  resolution 
study  was  conducted  to  obtain  reliable  heating  predictions.  Species  mass  fraction  and  non-dimensional  temper¬ 
ature  pro  les  along  the  stagnation  line  are  shown  in  Fig.  13  for  100  and  200  equidistant  point  grids,  for  both  the 
McKenzie  [Fig.  13(a)]  and  Park  [Fig.  13(b)]  bulk  chemistry  models.  The  solutions  for  both  mesh  re  nements 
are  seen  to  be  identical,  showing  the  grid  independence  of  the  results.  The  analysis  of  Fig.  13  reveals  a  non 
negligible  in  uence  of  the  bulk  chemistry  model.  The  main  difference  in  the  predictions  obtained  with  the 
two  sets  of  nite  rate  chemical  processes  lies  in  the  evolution  of  atomic  carbon  along  the  stagnation  line.  The 
McKenzie  model  does  not  enhance  considerably  C  recombination,  while  for  the  Park  model  C  is  involved  in 
the  formation  of  CO2  starting  from  80%  of  the  stagnation  line  and  it  is  completely  depleted  when  we  reach  the 
40%  of  the  stagnation  line.  As  a  consequence,  in  Fig.  13(a)  we  observe  that,  even  if  CO  disappears  completely, 
a  lower  amount  of  CO2  is  present  at  the  wall  with  respect  to  the  computations  obtained  with  the  Park  model. 
In  addition,  since  less  carbon  is  available  for  recombination,  a  higher  amount  of  O2  is  observed  at  the  wall  for 
the  prediction  obtained  with  the  McKenzie  model.  The  difference  in  the  prediction  of  atomic  carbon  chemistry 
will  be  more  evident  for  the  high  enthalpy  case  and  some  of  the  previous  observations  will  be  helpful  for  the 
de  nition  of  the  set  of  wall  reactions. 

We  now  turn  our  attention  to  the  analysis  of  the  diffusion  demixing  phenomenon  considering  three  different 
wall  boundary  conditions  : 
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Figure  14:  Mass  fractions  and  non-dimensional  temperature  for  three  different  boundary  conditions  (Park  bulk 
chemistry  model)  [ hs=  I 5.3MJ/kg.  FH  =300K|. 


•  fully  catalytic  wall  (present  wall  chemistry  model  -  FCW); 

•  equilibrium  wall  with  the  element  fraction  of  the  outer  edge  (LEW  -  A'); 

•  equilibrium  wall  with  the  element  fraction  corresponding  the  fully  catalytic  case 
(LEW  -  XcFC). 


The  wall  heat  ux  corresponding  to  these  cases  is  summarized  in  Tab.  6. 


Table  6:  Wall  heat  ux  for  several  boundary  conditions  (Park  chemistry  model)  [hg=  1 5.3MJ/kg,  rw-300K]. 


Boundary  Condition 

qw\MW/m2\ 

FCW 

0.856 

LEW  -  XcFC 

0.857 

LEW  -  A/ 

0.918 

Assuming  a  local  equilibrium  wall  with  the  outer  edge  elemental  fraction  results  in  an  overestimation  of  the 
wall  heat  ux,  which  is  explained  as  follows.  The  wall  heat  ux  is  primarily  controlled  by  the  recombination 
processes  taking  place  at  the  wall.  Due  to  diffusion  demixing  in  the  boundary  layer,  there  is  more  oxygen 
and  less  carbon  at  the  wall  then  at  the  outer  edge  (see  Fig.  14  which  shows  the  mass  fraction  pro  les  for  the 
3  boundary  conditions).  As  a  result,  fewer  CO2  molecules  can  be  formed,  the  excess  oxygen  forming  CL 
molecules.  Hence,  less  heat  is  released  than  if  C  and  O  elements  were  in  the  1/3  -  2/3  proportion. 

These  results  point  out  the  importance  of  diffusion  demixing  effects,  in  particular  regarding  their  in  uence  on 
wall  heat  ux,  and  clearly  establishes  that  the  local  equilibrium  boundary  condition  based  on  outer  edge  (free 
stream)  elemental  fractions,  widely  used  in  the  framework  of  hypersonic  o  ws  simulations,  is  incorrect  [61,  60]. 
At  this  point,  it  is  interesting  to  comment  on  the  relation  between  elemental  diffusion  and  surface  chemistry. 
Indeed,  some  approximate  models  used  to  describe  stagnation  line  o  ws  are  based  on  the  constant  elemental 
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fraction  assumption  as  for  example  in  the  y-codc  developed  at  IPM  [36].  Under  chemical  non  equilibrium 
conditions,  using  the  2-surface  reactions  model,  both  CU  and  CO2  should  be  present  at  the  wall  under  fully 
catalytic  conditions.  On  the  other  hand,  if  constant  elemental  fraction  is  assumed  through  the  boundary  layer, 
equal  to  the  outer  edge  one  (1/3-2/3),  the  co-presence  of  these  two  species  is  not  possible.  Indeed,  as  shown  in 
Ref.  [66],  the  results  of  the  y-codc  operating  with  CO2  mixtures  under  non  equilibrium  and  for  a  fully  catalytic 
surface,  based  on  the  2-surface  reactions  scheme,  show  the  presence  of  only  CO2  on  the  surface  for  cold  wall 
conditions.  This  is  clearly  in  contradiction  with  the  equal  effective  recombination  probability  which  is  one  of 
the  major  assumption  of  the  IPM  methodology.  We  wish  to  stress  that  the  model  presented  in  Sec.  5.1.1  is  free 
of  these  inconsistency. 

Another  illustration  of  the  importance  of  the  demixing  phenomenon  is  provided  by  comparing  Scott’s  and  the 
present  formulation  for  a  fully  catalytic  wall.  Computed  elemental  fraction  pro  les  are  presented  in  Fig.  15, 
always  for  the  low  enthalpy  case.  Elemental  fractions  of  carbon  (Xc)  and  oxygen  (X°)  de  ned  in  Eq.  (5),  are 
related  to  molar  elemental  concentration  (£,■)  as  follows  : 


X° 


£ o 

(rc  +  £ o 


and  Xc 


£c  +  £ o 


where 

£c  =  (xco2  +  xco  +  xc)  -j—  and  t;o  =  (2 xco2  +  xco  +  2 xo2  +  xo) 

Ru  being  the  universal  gas  constant.  From  the  analysis  of  Fig.  15(a)  we  observe  that  for  Scott’s  formulation 
X°  >  X°’6  throughout  the  boundary  layer,  which  is  in  contradiction  with  the  elemental  balance  [59].  Indeed, 
in  general  only  convective  ux  es  are  signi  cant  in  the  elemental  balance  (the  wall  diffusion  ux  is  zero  and  the 
edge  diffusion  ux  is  negligible)  so  that  the  elemental  balance  [59]  reduces  to 


tjsVs  + 


r,dUS 

dx 


5 


TJmax 

f  dR 
J  p 
0 


Fin) 

p 


dfj  -  0. 


(22) 


Hence,  the  elemental  concentration  (gj)  distribution  must  be  such  that  some  weighted  average  of  it  equals  the 
outer  edge  concentration,  and  similarly  for  the  elemental  fraction  (X7),  which  is  not  veri  ed  when  X°  >  X° -6 
throughout  the  boundary  layer.  In  contrast,  the  elemental  fraction  distributions  obtained  with  the  present  wall 
chemistry  model  are  compatible  with  this  condition  [Fig.  15(b)]. 

In  fact,  the  elemental  fraction  distribution  obtained  with  Scott’s  wall  chemistry  model  violates  this  condition 
because  this  model  does  not  respect  elemental  conservation  at  the  wall.  For  the  present  problem,  it  provides  a 
source  of  O  and  a  sink  of  C  at  the  wall.  This  was  veri  ed  by  computing  the  elemental  diffusion  ux  es  at  the 
wall  and  observing  that 

p  p  P  p  P  p 

q  .  CP2w  +  J  COw  +  JCw  <  q  q  .  2  C°2W  +  COw  +  2  °2W  +  Ow  y  Q 

Xfco2  Mco  Me  Mco2  Mco  M0l  Mo 

We  now  analyze  the  calculated  heat  ux  maps  for  the  present  o  w  conditions.  Heat  ux  maps  consist  of  curves 
of  wall  heat  ux  iqw)  as  a  function  of  wall  temperature  between  300  K  and  2100  K  for  various  values  of  the 
recombination  probability  yw.  Four  different  maps  have  been  obtained  for  the  present  test  case,  corresponding 
to  all  possible  bulk  chemistry  (Park/McKenzie)  and  wall  chemistry  (Scott/present)  model  combinations.  They 
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(a)  Scott  wall  chemistry  model. 


(b )  Present  wall  chemistry  model. 


Figure  15:  Element  fractions  computed  with  the  Scott’s  and  the  present  catalycity  models  (Park  bulk  chemistry 
model)  [0:  wall;  1:5]. 


(a)  McKenzie  bulk  chemistry  model. 


(b)  Park  bulk  chemistry  model. 


Figure  16:  Heat  ux  maps  for  the  lower  enthalpy  case  [/t<s=15.3MJ/kg]. 


are  displayed  in  Fig.  16. 

The  in  uence  of  the  bulk  chemistry  model  is  clearly  visible.  Park’s  model  reaction  rates  being  much  larger 
than  McKenzie’s,  the  wall  heat  ux  is  signi  cantly  larger  for  low  wall  temperatures  and  low  recombination 
probabilities.  The  in  uence  of  the  wall  reaction  model  is  also  important  for  large  recombination  probabilities. 
For  high  values  of  catalycity,  the  heat  ux  is  seen  to  be  signi  cantly  higher  with  the  present  model,  which  is 
consistent  with  the  fact  that,  because  of  the  existence  of  a  sink  of  C  atoms  at  the  wall  in  Scott’s  model,  there 
are  fewer  C  atoms  available  for  recombination. 
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5.1.3  High  enthalpy  case 


Attention  is  now  turned  to  the  high  enthalpy  case.  Because  of  the  higher  enthalpy  level,  the  mass  fraction  of 
C  atom  at  the  outer  edge  is  much  higher,  so  that  the  third  wall  reaction,  C  +  2(0  -  S)  —>  CO 2.  may  become 
important.  For  this  reason,  we  rst  analyze  the  in  uence  of  the  wall  reactions  set  using  Scott’s  formulation 
for  both  bulk  chemistry  models.  Species  mass  fractions  and  non-dimensional  temperature  pro  les  are  shown 
in  Fig.  17  for  a  fully  catalytic  wall.  Whereas  for  Park’s  model  there  is  no  in  uence  of  the  wall  reactions  set 
because  the  carbon  atom  entirely  recombines  inside  the  boundary  layer  (bulk  chemistry),  on  the  contrary  the 
wall  reactions  set  has  a  large  in  uence  when  using  McKenzie’s  model,  for  which  bulk  reactions  are  much 
slower. 


(a)  McKenzie  bulk  chemistry  model.  (b)  Park  bulk  chemistry  model. 


Figure  17:  Mass  fractions  and  non-dimensional  temperature  pro  le  [ /i,»,-39MJ/kg,  7\V=300K][0:  wall;  1:  5]. 


In  particular,  if  the  third  wall  reaction  is  ignored,  there  remains  an  important  amount  of  atomic  carbon  at  the 
wall.  As  a  result,  the  wall  heat  ux  is  considerably  smaller  (see  Tab.  7  in  which  heat  ux  es  are  listed  for  all 
bulk  chemistry/wall  reaction  set  combinations). 


Table  7:  Wall  heat  ux  for  several  sets  of  surface  reactions  [/t£=39MJ/kg,  7’H,=300K  ]. 


Surface 

Stagnation  Line 

qw[MW/m2] 

Scott  2WR 

McKenzie 

1.20 

Scott  3WR 

McKenzie 

2.14 

Scott  3WR 

Park 

2.23 

Scott  2WR 

Park 

2.24 

The  present  formulation  for  the  3  wall  reactions  model  reads  as  follows 


7o  +  Vo  +  Vo  =  vo,  y2co  =  vco,  yl  =  re.  y20Ml0  =  y2coH:o . 


y3cMlc  = 


v*(M) 

2 


(23) 
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where  yco  —  7o  =  7c  =  Jw  Eqs.  (23)  represent  a  system  of  5  equations  in  5  unknowns  ylQ,  y2(),  y3Q,  y2co,  and 
y3,.  For  the  present  o  w  conditions,  numerical  experiments  have  shown  that  the  solution  does  not  satisfy  the 
constraint  y3Q  <yw,  which  implies  that  the  hypothesis  of  equal  recombination  probabilities  for  all  species  does 
not  hold  [7 1].  The  analysis  of  a  suitable  assumption  is  however  left  as  an  open  subject  for  a  further  study. 

We  therefore  switch  our  attention  to  the  Park  bulk  chemistry/2  wall  reaction  model  for  the  present  high  enthalpy 
case. 


(a)  Mass  fractions  and  non-dimensional  temperature 
pro  le.  (r„,=300K)  (0:  wall; 


Figure  18:  Stagnation  line  solution  for  the  higher  enthalpy  case  (Park  bulk  chemistry  model)  [  /t^-39MJ/kg]. 

Species  mass  fraction  and  non-dimensional  temperature  pro  les  for  a  fully  catalytic  wall  computed  using  the 
present  wall  chemistry  formulation  are  shown  in  Fig.  18(a).  This  is  to  be  compared  with  Fig.  17(b)  where 
Scott’s  formulation  was  used.  The  same  conclusion  as  for  the  low  enthalpy  case  is  obtained,  i.e.  for  high  yw 
Scott’s  model  acts  as  a  sink  of  carbon  atoms  at  the  wall.  As  a  result,  the  wall  heat  ux  is  smaller,  as  seen  in 
Fig.  18(b)  which  shows  the  heat  ux  map  for  the  high  enthalpy  case. 


5.1.4  In  uence  of  outer  edge  elemental  fractions  on  heat  ux  maps 


As  discussed  in  the  introduction,  the  methodology  developed  at  the  IPM  and  currently  used  at  VKI  to  estimate 
the  catalycity  properties  of  TPS  materials,  is  based  upon  Navier-Stokes  simulations  of  an  ICP  facility  supposed 
to  be  under  FTE  conditions  determined  assuming  constant  elemental  fractions. 

The  puipose  of  the  Navier-Stokes  plasma  o  w  simulation  is  to  provide  the  FTE  outer  edge  inputs  for  the  nite 
thickness  non  equilibrium  stagnation  line  o  w  computations,  earned  out  to  compute  a  heat  ux  map.  Among 
the  quantities  needed,  the  composition  is  implicitly  given  as  a  function  of  the  outer  edge  pressure  and  rebuilt- 
enthalpy,  assuming  FTE  and  using  the  torch  inlet  elemental  fractions.  In  general,  thanks  to  elements  diffusion, 
the  elemental  fractions  vary  within  the  ICP  facility  [76],  as  shown  in  Sec.  4,  leading  to  an  outer  edge  elemental 
composition  different  from  the  torch  inlet  one.  If  this  difference  has  an  in  uence  on  the  heat  ux  map,  by  shift¬ 
ing  or  stretching  the  iso-catalycity  lines,  the  same  in  uence  will  be  re  ected  to  the  estimation  of  TPS  material 
catalycity. 
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ii/ii 


max 


(c)  Species  mass  fraction  :  CO>  mixture  ( T„ .  =  300X,  yw  = 
1)  (0:w;l:<?) 


(d)  Species  mass  fraction  :  air  mixture  ( T„  =  300  A'.  yw  = 
1)  (0:w;l:<5) 


Figure  19:  In  uence  of  outer  edge  elemental  fraction  on  the  stagnation  line  solution  [The  outer  edge  conditions 
are  speci  ed  in  Ref.  [58]]. 


One  way  to  easily  estimate  this  in  uence  is  to  use  the  results  of  the  constant  elemental  fraction  ICP  simulation 
and  vary  only  the  outer  edge  elemental  fractions,  used  to  determine  the  outer  edge  composition  for  the  given 
enthalpy  and  pressure.  In  Sec.  4  we  found  an  oxygen  elemental  fraction  lower  than  the  inlet  value  (X°  =  0.21) 
in  the  jet  near  the  axis  for  air  plasmas.  On  the  other  hand,  for  carbon  dioxide  o  ws  the  opposite  behavior  was 
observed,  i.e.  the  oxygen  elemental  fraction  reached  values  higher  than  at  the  torch  inlet  (X°  =  2/3).  Following 
this  result,  we  have  perturbed  the  oxygen  outer  edge  elemental  fraction  by  +10%  and  -10%  of  the  reference 
inlet  value,  for  carbon  dioxide  and  air  plasmas  respectively. 

The  results  of  this  analysis  are  presented  in  Figs.  19(a)  and  19(b),  where  are  shown  the  heat  ux  maps  corre¬ 
sponding  to  the  conditions  de  ned  in  Ref.  [58],  for  both  air  and  carbon  dioxide.  The  results  show  an  important 
in  uence  on  the  stagnation  point  heat  ux  for  carbon  dioxide  mixtures.  The  reason  of  these  differences  can  be 
explained  from  the  analysis  of  species  concentration  showed  in  Fig.  19(c)  for  a  fully  catalytic  cold  wall.  Indeed, 
the  increase  of  oxygen  outer  edge  fraction  reduces  the  amount  of  carbon  available  for  CO2  recombination  at  the 
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surface,  this  leading  to  a  lower  heat  ux.  This  is  more  evident  for  high  values  of  recombination  probabilities, 
where  almost  the  same  effect  on  wall  heat  ux  is  observed  in  the  whole  temperature  range  considered.  As  the 
recombination  probability  is  reduced  the  in  uence  of  the  outer  edge  elemental  fraction  decreases  as  well,  and  a 
higher  in  uence  is  observed  for  cold  wall  conditions,  where  recombination  is  enhanced  by  the  bulk  chemistry 
owing  to  the  low  temperature. 

As  far  as  air  mixtures  are  concerned,  the  effects  on  wall  heat  ux  arc  essentially  negligible.  Indeed  the  lines  of 
Fig.  19(b)  are  superposed  for  all  values  of  catalycity  in  the  temperature  range  investigated.  The  quite  low  in  u- 
ence  on  heat  ux  can  be  explained  by  the  analysis  of  the  species  mass  fraction  pro  les  of  Fig.  19(d).  Indeed  the 
perturbation  of  the  elements  fraction  shifts  slightly  the  concentration  of  N2,  N,  and  O  at  the  outer  edge  without 
altering  too  much  the  stagnation  line  chemistry. 

The  analysis  carried  out  in  this  section  shows  therefore  that,  especially  for  Mars  entry  applications,  the  outer 
edge  elemental  fractions  play  an  important  role  in  the  determination  of  the  heat  ux  map  and  therefore  in  the 
extrapolation  of  the  TPS  material  catalycity. 


6  Conclusions 


In  this  section  the  main  achievements  of  the  present  research  are  recalled  following  the  structure  of  this 
manuscript  with  the  aim  to  summarize  the  contributions  brought  to  the  modeling  of  high  enthalpy  gases  in 
the  framework  of  a  macroscopic  methodology  used  ti  estimate  the  catalytic  properties  of  thermal  protection 
materials. 


6.1  Theory  of  reacting  mixtures  under  Local  Thermodynamic  Equilibrium 

A  series  of  new  elemental  transport  coefficients  have  been  presented  for  mixtures  under  thermochemical  equi¬ 
librium  conditions.  They  allow  to  reduce  the  equations  of  chemically  reacting  o  ws  in  LTE  to  an  elegant 
system  consisting  of  the  conventional  Navier-Stokes  equations  (mass,  momentum,  energy)  complemented  by 
an  advection-diffusion  equation  for  the  mass  fraction  of  each  chemical  element  in  the  mixture.  The  obtained 
formalism  is  in  closed  form  in  the  sense  that  diffusive  ux  es  are  directly  expressed  in  terms  of  gradients  of 
the  solution  unknowns,  unlike  other  formulations  in  which  these  ux  es  are  obtained  in  an  implicit  manner,  by 
solving  the  full  system  of  Stefan-Maxwell  equations.  In  the  case  of  mixtures  containing  neutral  components, 
assuming  local  neutrality  and  absence  of  current,  we  have  introduced  the  following  new  transport  properties  : 
the  elemental  multicomponent  ( Deq )  and  thermal  demixing  (Dj)  diffusion  coefficients,  the  thermal  demixing 
conductivity  (Ad),  the  elemental  heat  transfer  coefficients  (AeEI).  We  refer  the  interested  reader  to  Ref.  [57]  for 
further  details  about  the  extension  of  this  formalism  composed  of  both  neutral  ad  charged  components. 
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6.2  Elemental  diffusion  along  the  stagnation  line 

To  verify  the  importance  of  elemental  diffusion  phenomena  and  to  assess  the  performances  of  the  theory  pre¬ 
sented  in  Sec.  2,  we  focused  on  a  8-species  CO2/N2  mixture  containing  three  elements  (C,  N,  O).  This  analysis 
was  based  on  the  use  of  the  transport  properties  previously  presented,  which  required  the  computation  of  nine 
elemental  multicomponent  and  three  thermal  demixing  coefficients,  the  demixing  thermal  conductivity  and 
three  elemental  heat  transfer  coefficients,  in  addition  to  the  ‘standard’  LTE  transport  properties. 

The  results  of  this  study  represent  the  rst  attempt  to  compute  thermal  demixing  and  multicomponent  diffusion 
coefficients  as  well  as  elemental  heat  transfer  coefficients  and  thermal  demixing  conductivity  for  a  mixture  of 
heteronuclear  components  containing  three  elements.  This  proves  the  generality  of  the  proposed  formulation 
and  clearly  represents  an  improvement  with  respect  to  previous  formulations  available  in  the  literature.  In  ad¬ 
dition,  the  results  obtained  with  the  proposed  closed  form  of  the  equations  coincide  with  those  obtained  using 
an  implicit  formulation,  indicating  the  correctness  of  the  new  formulation. 

Finally,  after  a  critical  analysis  of  the  results,  we  noticed  how  complex  elemental  concentration  and  heat  ux 
patterns  can  be  easily  explained  by  examining  graphs  of  LTE  transport  coefficients  that  arise  from  the  new 
proposed  LTE  formulation. 


6.3  Elemental  diffusion  in  inductive  coupled  plasmas 

In  Sec.  4,  we  extended  our  analysis  of  elemental  diffusion  considering  the  investigation  of  air  and  carbon 
dioxide  inductively  coupled  plasma  o  ws  for  application  to  Earth  and  Mars  entry.  To  verify  the  validity  of 
one  of  the  major  hypothesis  of  the  IPM  methodology,  i.e.  the  existence  of  equilibrium  conditions  within  ICP 
facilities,  we  considered  both  the  equilibrium  and  the  non  equilibrium  form  of  the  governing  equations.  The 
numerical  solution  of  these  equations  led  to  a  detailed  analysis  of  the  o  w  behavior  within  the  torch  for  air  and 
CO2  mixtures  and  then  of  the  jet  o  wing  into  the  test  chamber  for  air  mixtures. 

Firstly,  we  presented  original  carbon  dioxide  plasma  o  w  simulations  under  equilibrium  considering  both 
constant  and  variable  elemental  fraction  conditions.  A  critical  analysis  of  these  results  shows  that  for  an  inlet 
carbon  volumetric  fraction  of  ~  33%  : 


•  the  solution,  under  chemical  equilibrium,  of  an  additional  set  of  elemental  continuity  equations  for  oxy¬ 
gen  and  carbon  shows  the  presence  of  important  variations  in  the  elemental  composition. 

•  The  amount  of  demixing  observed  for  the  operating  conditions  considered  is  characterized  by  a  maximum 
increase  to  ~  38%  and  a  decrease  to  ~  28%  of  the  carbon  elemental  fraction. 

•  The  elfects  of  these  variations  in  the  elemental  fraction  on  the  species  evolution  is  much  higher  than  for 
air  o  ws.  Indeed,  at  the  outlet  of  the  torch,  important  differences  are  observed  between  the  species  radial 
pro  les  obtained  under  equilibrium  with  constant  and  variable  elemental  fractions. 

•  The  in  uence  of  demixing  on  the  temperature  prediction  is  much  higher  for  CO2  o  ws  than  for  air.  At  the 
torch  exit,  neglecting  elemental  diffusion  a  maximum  relative  error  of  ~  40%  is  observed  for  r  -5.8  cm, 
while  on  the  axis  it  decreases  sharply  to  ~  2%. 
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Secondly,  we  turned  our  attention  to  the  analysis  of  air  plasma  o  ws  both  in  the  torch  and  in  the  chamber 
of  the  VKI-Plasmatron.  Two  operating  pressures  have  been  investigated  corresponding  to  5  kPa  and  30  kPa. 
In  addition,  the  three  chemical  regimes  used  to  analyze  air  plasmas  within  the  torch  have  been  considered 
and  for  the  description  of  chemical  non  equilibrium,  two  nite  rate  chemistry  models  have  been  selected  as 
representative  of  Earth  atmosphere  chemistry. 

The  analysis  of  these  results  leads  us  to  the  following  conclusions  : 


•  for  the  low  pressure  case  (5  kPa),  we  observe  a  certain  in  uence  of  the  nite  rate  chemistry  model  on  the 
prediction  of  the  o  w  behavior  within  the  torch,  while,  as  we  move  towards  the  test  chamber,  the  o  w 
predictions  rapidly  match  leading  to  almost  negligible  discrepancies  ~  30  cm  after  the  torch. 

•  Non  equilibrium  effects  are  observed  within  the  torch  for  the  5  kPa  operating  pressure.  In  the  jet,  the  non 
equilibrium  results  are  rather  close  to  the  equilibrium  ones,  provided  that  elemental  demixing  is  taken 
into  account. 

•  The  consequence  of  increasing  the  operating  pressure  is  to  sharply  decrease  the  differences  between  the 
equilibrium  and  non  equilibrium  predictions.  These  differences  are  very  small  for  the  low  pressure  case 
(5  kPa)  and  they  become  negligible  for  the  high  pressure  case  (30  kPa),  especially  in  the  jet. 


The  major  result  of  this  analysis  consists  in  the  fact  that,  for  the  purpose  of  TPM  testing,  the  costly  sim¬ 
ulation  of  chemical  non  equilibrium  ICP  o  ws,  both  in  the  torch  and  in  the  test  chamber,  can  be  avoided  for 
sufficiently  high  pressures.  Indeed  the  same  accuracy  can  be  reached  relying  on  the  numerical  simulation  of  a 
more  compact  thermochemical  equilibrium  formalism  based  on  the  solution  of  an  additional  set  of  elemental 
advection  diffusion  equations,  free  of  stiff  source  terms.  In  addition,  because  of  the  important  elemental  frac¬ 
tion  variations  observed  as  a  result  of  both  air  and  carbon  dioxide  plasma  o  w  simulation,  it  appears  desirable 
to  add  the  elemental  fractions  at  the  stagnation  line  outer  edge  to  the  list  of  non  dimensional  parameters  de  ned 
in  the  framework  of  the  IPM  methodology. 


6.4  Practical  application  to  TPM  testing 

Following  the  results  obtained  in  the  previous  sections,  we  focused  our  interest  on  the  analysis  of  stagnation 
line  o  ws  for  application  to  TPS  testing.  To  this  end,  both  air  and  carbon  dioxide  mixtures  have  been  analyzed 
considering  only  neutral  components. 

In  our  simulations,  two  models  (Park  [52]  and  McKenzie  [44])  for  the  description  of  nite  rate  chemistry  for 
a  carbon  dioxide  mixtures  of  neutral  species  have  been  tested  and  important  differences  in  the  results  have  been 
observed.  Then,  our  attention  has  been  attracted  by  the  heterogeneous  recombination  processes  taking  place 
on  the  surface  of  TPM  samples.  During  our  investigation,  some  unexpected  a  ws  of  the  Scott  [69]  formulation 
have  been  discovered.  To  overcome  them,  an  alternative  model  has  been  proposed  and  successfully  tested. 

The  main  outcome  of  the  analysis  of  the  results  obtained  with  this  new  gas/surface  interaction  model  reads  as 
follows  : 
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•  a  deeper  description  of  the  theoretical  relation  between  a  fully  catalytic  surface  and  the  related  equilib¬ 
rium  conditions  is  presented; 

•  we  have  shown  that  the  constant  elemental  fraction  assumption  may  contradict  the  hypothesis  of  equal 
effective  recombination  probability; 

•  we  have  shown  that  a  non  equilibrium  model  of  a  carbon  dioxide  stagnation  line  o  w  based  on  the 
constant  elemental  fraction  assumption  is  in  contrast  with  the  equal  recombination  probability  for  CO 
and  O  involved  in  the  formation  of  CO2  and  O2  at  the  wall; 

•  we  have  shown  that  the  equal  recombination  probability  hypothesis  may  not  be  compatible  with  the  local 
wall  conditions; 

•  we  have  veri  ed  the  in  uence  of  the  outer  edge  elemental  fractions  on  the  stagnation  point  heat  ux.  This 
in  uence,  rather  important  for  carbon  dioxide  mixtures,  con  rms  the  necessity  to  add  the  them  to  the  list 
of  parameters  of  the  methodology  developed  at  the  IPM  suited  for  the  estimation  of  catalytic  properties 
of  TPS  materials,  especially  for  Mars  entry  applications. 
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